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ELEMENTS  OF  FREE-STREAMLINE  THEORY 

by 

•Joiix  S.  McXowx* 


J X TKODU CTOKY  REMARKS 


Many  plane  irrotational  flow  patterns  which  include  the  for- 
mation or  deflection  of  free  jets  can  be  analyzed  completely  by  a 
method  known  as  the  Helmholtz-Kirchhoff  theory  of  free  stream- 
lines. The  results  obtained  may  be  significant  not  only  for  the 
corresponding  two-dimensional  flows  of  real  incompressible  fluids 
but,  in  some  instances,  for  their  three-dimensional  counterparts 
as  well.  Numerous  worthwhile  calculations  remain  to  be  made, 
perhaps  because  the  theoretician  and  the  hydraulieian  have  not 
yet  pooled  their  interests  for  effective  exploitation  of  this  powerful 
tool.  The  analyses  consist  of  the  definition  of  successive  conformal 
transformations  involving  a hodograph,  or  velocity  plane,  and  the 
application  of  the  Schwarz- Oh ristoffel  transformation.  A number 
of  such  calculations  are  presented  in  the  following  articles. 

Impetus  for  the  preparation  of  this  bulletin  came  from  the 
writer’s  earlier  contact  with  Professor  AY.  E.  Brooke  (Chairman 
Emeritus  of  the  Mathematics  and  Mechanics  Department,  Univer- 
sity of  Minnesota),  and  from  suggestions  by  Professor  Kravtchenko 
(Faculte  d s Sciences,  Universite  dc  Grenobie)  and  M.  DuPort 
(Engineer  at  the  Etablisscments  Xcyrpie,  also  in  Grenoble).  The 
second  paper  was  prepared  by  M.  Gerber  and  the  writer,  at 
Kravtchenko's  suggestion,  in  order  to  obtain  a comparison  of  the 
results  of  theory  with  those  from  experiments  with  cavitating 
flows.  The  remaining  papers  were  prepared  subsequently  by 
second-  and  third-year  graduate  students  in  the  department  of 
Mechanics  and  Hydraulics  at  the  State  University  of  Iowa.  This 
cooperative  endeavor  was  undertaken  as  an  experiment  in  com- 
bining teaching  with  the  conduct  of  original  investigation.  These 
selecte  i students,  many  of  whom  are  also  staff  members  of  the 
Iowa  institute  of  Hydraulic  Research,  began  with  a literature 
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survey  and  group  review  of  the  essential  tlieory.  and  then  eum- 
pleted  the  detailed  analyses  and  calculations  of  the  various  prob- 
lems selected.  Formal  oral  reports  were  coordinated  and  presented 
to  an  audience  composed  of  faculty  members  and  other  graduate 
students,  end  the  pedagogical  experiment  was  completed  by  the 
preparation  of  the  remainder  of  the  papers  in  this  bulletin.  Mr. 
Carstens  of  the  Department  of  Engineering  English  provided 
valuable  criticism  and  assistance  in  their  preparation.  Financial 
assistance  was  provided  by  the  Office  of  Naval  Research  (NSonr- 
f>Of)  TO  •rD  for  the  conduct  of  the  experiments  and  for  the  editing 
and  publication  of  this  Bulletin. 

Historical  Dkvki.oujik.nt 

During  the  period  ot  over  eighty  years  which  has  elapsed  since 
the  original  developments  of  Helmholtz  [1]*  and  Kirehhoff  [2],  the 
numerous  contributions  to  Ibis  intriguing  phase  of  analytic  function 
theory  have  been  made  primarily  by  mathematicians,  with  results 
which  sometimes  appear  bizarre  from  the  practical  point  of  view. 
Aside  from  a few  classic  cases  which  were  solved  before  1900  [1. 
2.  ;V|.  the  kinds  of  flows  analyzed  have  usually  been  significant  more 


Fin.  1.  Com r.Niusnx  oe  Kvtkam  k Coniutjons 
(a)  No  Skpahation;  (b)  Fki:i:  Stiu:.\mmnk 


for  the  mathematical  method  than  for  the  result.  The  mathematician 
has  seldom  been  interested  in  the  applications  of  tlm  subject,  where- 
as t lie  iiydrauiiciau  with  much  to  gain  has  equally  rarely  acquaint- 
ed himself  with  the  tlieory  sufficiently  well  oven  to  make  sugges- 
tion* to  the  theoretician.  Tints  with  only  occasional  exceptions,  such 
as  the  works  of  von  Mises  [4]  and,  much  more  recently,  Birkhoff 
et  al  }">]  and  Plcsset  and  Shaffer  [61.  little  of  value  to  the  prac- 

‘■Kcli't ■ i' a ces  for  all  of  the  pavers  are  collected  at  the  end  of  this  Bulletin. 
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titionor  has  been  presented  in  the  Humorous  works  on  this  subject 
since  publication  of  the  earliest  papers. 

Helmholtz  devised  tnc-  concept  of  the  free  streamline  in  a search 
for  a more  realistic  solution  to  the  problem  of  efflux  from  (or  flow 
into)  a two-dimensional  pipe  or  channel.  The  only  pr:or  classical 
theory  was  for  flow  which  doubled  back  on  the  boundary  in  such 
a way  as  to  give  infinite  velocities  and  negatively  infinite  pressures 
at.  the  end  of  the  channel  (point  A in  Fig.  la).  His  fundamental 
contribution  was  the  concept  of  a free  boundary  which  was  defined 
in  the  kinematic  rather  than  in  the  geometric  sense  (Fig.  lb).  He 


Fig.  ‘1.  Fi.ow  Patterns  for  (a)  Slot,  and  (b)  Flat  Plati: 

reasoned  logically  that  the  bounding  streamline  would  separate 
from  the  solid  boundary,  and  that  the  free  streamline  thus  formed 
could  be  characterized  by  a constant  pressure  and,  hence,  by  a 
constant  velocity.  He  visualized  a quiescent  wake  of  constant  pres- 
sure and  a velocity  discontinuity  at  the  free  streamline.  "With  the 
aid  of  this  concept,  he  proceeded  to  derive  the  equations  for  the 
free  jet  issuing  from  a Bovda  tube.  Kirehhoff  shortly  afterwards 
devised  functions  for  the  description  of  the  efflux  from  an  opening 
in  an  infinite  reservoir  having  plane  boundaries,  and  of  the  effect 
of  a plate  on  a field  of  flow  which  is  otherwise  uniform  and  extends 
to  infinity  in  all  directions  (Fig.  2).  Rayleigh  1 3 ] systematized 
these  results  somewhat  and  added  his  analysis  of  intersecting  jets. 
The  method  employed  by  these  pioneers  was  indirect  in  that  it 
consisted  of  trying  various  functions  and  adjusting  constants  until 
some  useful  result  was  obtained. 

Development  of  a direct  method  of  solving  any  of  a large  class 
of  free-streamline  flows  awaited  definition  of  the  hodograph  planes 
and  certain  complementary  transformations.  The  ordinary  hodo- 
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graph  or  vdoeiiy  piano  was  first  used  by  Miolieil  |_T J in  conjunc- 
tion with  the  well  Known  transformation 


Planck  IS]  in  1894  suggested  the  logarithmic  hodograph  which, 
together  with  the  important  transformation  theorem  developed  by 
Christoffel  [9]  in  1867  and  independently  by  Schwarz  [10]  in 
I860,  made  possible  the  direct  resolution  oi  flow  patterns  for  a 
wide  variety  of  houmlary  configurations.  The  primary  restrictions 
on  the  method,  aside  from  complexity  of  detail,  were  the  general 
considerations  that  the  motions  he  two-dimensional,  incompressible, 
and  unaffected  by  gravity  or  viscosity,  and  that  the  solid  boundaries 
be  composed  of  rectilinear  segments. 

In  the  interval  between  1890  and  1920  a number  of  contributions 
were  made  of  flow  patterns  for  a variety  of  boundary  conditions. 
AVitli  the  exception  of  the  very  few  papers  on  the  simplest  and 
most  fundamental  problems  and  of  the  paper  by  von  Mises  [4], 
already  mentioned,  few  calculated  results  were  presented,  the 
authors  being  content  to  give  general  results  or  simply  to  indicate 
their  existence. 

Not  until  the  past  few  years  lias  there  been  any  renewal  of 
interest  in  the  extension  and  application  of  this  type  of  analysis 
to  obtain  useful  results  in  current  investigations.  Recently,  Birk- 
hoff,  Plesset,  and  others  have  published  studies  of  boundary  effects 
in  cavity  flows,  and  Salzman  [11]  and  McNown  and  Hsu  [12] 
have  published  analyses  of  lateral  efflux  through  an  opening  in 
one  side  of  a two-dimensional  channel.  Comparisons  of  the  results 
of  theory  and  of  experiment  have  indicated  a correspondence  which 
is  usually  dose  ami  sometimes  astonishingly  so. 


Thf.  Hodourami  Planes 


in  the  derivation  of  functional  relationships  which  make  pos- 
sible the  mathematical  description  of  a pattern  of  flow,  definition 
of  a hodograph  or  velocity  plane  is  an  essential  step.  Two  types 
are  used  which  differ  in  detail  but  are  similar  in  concept  in  that 
they  make  possible  the  comparatively  simple  representation  of  the 
complicated  hounding  streamlines.  These  two  planes,  on  which  the 
How  patterns  are  transformed  conformally  in  accordance  with  the 
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precepts  of  complex-function  funny,  can  host  bo  described  through 
presentation  in  connection  with  a specific  example. 

If  ic  — <p  -f  i if  represents  the  complex  potential  and  2 = ./•  -f  ill 
the  variable  in  the  plane  of  actual  flow,  the  equation 

v = f(z) 

indicates  tiie  relationship  sought.  A fundamental  concept  of  hydro- 
dynamics is  that 

dw 

— = — u + tr  = -or'9 
dz 


in  which  u-  and  v are  the  x-  and  //-components  of  the  velocity  at 
any  point , q is  the  magnitude  of  the  resultant  velocity,  and  0 is 
its  inclination  measured  from  the  g-axis,  as  shown  in  Fig.  3.  Tints, 


q = + v'v~  v-  . 9 = tan-1  — 

u 

If,  to  simplify  the  terminology,  £ is  substituted  for  div/dz,  the  two 
kinds  of  hodographs  commonly  used  are  representations  of  £ and 


y 


x 

Fir..  3.  Di  rixrnox  Skktcu  roa  Viaoen  v. 

of  if  = In  ( — Yj  £),  the  second  being  also  expressible  in  the  form 

o=\n(-V~)=  -In  4-  +i9 

V v d w / Vj 

Straight  solid  hoimdaires  in  the  2-plane  transform  into  vndial 
lines  in  the  £-planc  (u/v  = constant),  and  into  straight  lines  par- 
allel to  the  real  axis  (9  = constant)  in  the  ff-plane.  Furthermore, 
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free  streamlines,  along  which  the  pressure  and  velocity  are  con- 
stant in  the  --plane,  become  circular  arcs  with  centers  at  the  origin 
( >r  — v"  = constant)  in  the  c-planc,  and  straight  lines  parallel  to 
the  imaginary  axis  ( q = constant)  in  the  fi-plane.  It  is  thus  evident 
that  the  orginnl  boundary  transforms  into  either  a circular  sector 
or  a rectangle. 

As  an  illustration  of  this  method,  the  four  planes.  ( z , Z.  Cl,  and 
w)  are  shown  in  Pig.  4a  for  the  unsym.metrieal  efflux  from  a two- 
dimensional  channel,  corresponding  points  having  been  assigned 
the  same  letter  in  each  plane.  For  mathematical  simplicity,  a ref- 
erence velocity  such  as  Vj  is  often  assigned  the  value  unity.  This 
has  not  been  done  herein  so  that  the  physical  significance  of  the 
variables  would  not  be  altered.  The  points  A and  D represent 
regions  far  to  the  left  and  far  to  the  right,  respectively,  at  which 
the  velocity  is  constant.  The  <p-  and  i/'-lines  in  the  w;-plane  arc 
vertical  and  horizontal  as  always,  and,  for  this  case  of  finite  quan- 


Fig.  4.  Conformai.  Mapping  Planks. 
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tity  of  flow,  the  boundary  of  !he  flow  in  the  w-planc  becomes 
a rectangle  of  fin ite  width  (0  5=  i p^Y0b)  and  infinite  length 
(—  cc  5S  <£  ^ c c).  The  coordinates  of  the  various  points  are  indi- 
cated in  Table  1 ; the  /-plane  is  discussed  in  the  following  section. 

Table  I 


z 

c 

V.  w 

t 

A 

— x -f-  i (0  to  b) 

—Vo 

— ln(Vo/F;)  x>  -f -i  (0  to  Vob)  f = 

-l(YL+YL\ 
A1?  ’V  / 

B 

i b 

0 

-f-  i (0  to  — 7t/2)  4>B~ViVob 

uO 

C 

ia 

-iV, 

— i(x/2)  4>c  ViVob 

1 

D 

undefined  - 

-V,€l 

3 —i/3  — °c-K(0  to  Vob) 

g = —cos  2/3 

E 

0 

-Vi 

0 0 

-1 

As  the  regions  at  A and  D become  points  in  the  £-  and  O-planes, 
the  concept  of  flow  from  a point  source  to  a point  sink  v, itliin  pre- 
scribed boundaries  is  useful  for  the  interpretation  of  these  inter- 
mediate planes.  That  is.  they  are  not  only  plots  of  the  velocity, 
but  transformed  flows  as  well;  intermediate  streamlines  are  indi- 
cated in  the  several  planes.  The  necessary  simplification  of  the 
boundary  has  been  accomplished  by  liie  representation  in  either 
the  £-  or  the  fi-plane.  The  links  between  the  several  planes  remain 
to  be  defined. 

Definition  of  Transformations 

In  order  to  utilize  the  planes  in  the  definition  of  the  essential 
fu  letional  relationship  between  v : and  z,  either  ^ or  fi  must  be 
rela  ed  to  iv  by  means  of  known  methods  through  still  another 
complex  variable  designated  as  (.  This  relationship  completes  the 
rather  complicated  functional  analysis  because,  as  has  already  been 
pointed  out,  either  £ or  A can  be  expressed  in  terms  of  dw/dz.  The 
only  remaining  difficulties  arise  in  the  explicit  evaluation  of  the 
various  functions  which  are  initially  known  only  as  differential 
equations 

Any  polygonal  boundary  in  one  complex  plane  can  be  trans- 
formed into  the  entire  real  axis  of  another  through  application  of 
the  Schwarz-Uhristoffel  theorem,  the  interior  of  the  polygon  being 
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transformed  into  either  the  upper  or  the  lower  half  of  the  second 
plane  [13].  Thus  the  outline  in  the  12-plnne  and  that  in  the  n-plane 
are  directly  transformable  into  the  real  axis  of  the  (-plane  as  shown 
in  Fig.  4h.  The  transforming  function  can  be  written  as  follows: 

o = \{  I -x-  \ 

" J ((-(,)*>/*  (t-k)*’1’ 

As  many  factors  ((  — tr.)  are  introduced  into  the  transformation  as 
there  are  vertices  of  the  polygonal  boundary  in  the  12-plane  (see 
Fig.  5).  The  value  (»  is  the  locus  of  the  transformed  point  on  the 
real  axis  in  the  (-plane,  and  /?„  is  the  exterior  angle  of  the  polygon 
at  that  point  in  the  Q-plane.  Three  of  the  (,,-values  can  be  assigned 
arbitrarily  and  the  remainder  must  be  given  parametric  values  to 
be  evaluated  in  terms  of  given  quantities,  along  with  M and  N, 
from  the  integrated  function.  Points  for  which  i„  is  infinity  need 
not  he  included  no  matter  what  the  value  of  /?„  . If  the  interior 
angle  exceeds  tt-,  the  exterior  angle,  as  well  as  the  corresponding 
exponent,  becomes  negative. 

In  defining  the  transformation  between  12  and  ( for  the  example 
of  Fig.  4,  the  points  B,  C,  and  E arc  arbitrarily  placed  in  the 
(-plane  at  oo,  -4-1,  and  — 1,  respectively.  The  unknown  or  para- 
metric values  of  the  abscissas  at  A and  D arc  assigned  the  values 
/ and  g.  The  vertices  of  the  polygon  are  at  the  points  B,  C,  and  E 
for  which  the  exterior  angles  are  tt,  tt/2,  and  tt/2.  Because  B is 
placed  at  oo,  the  only  terms  to  be  included  are  those  at  C and  E, 
and  the  transformation  becomes  a familiar  one, 

°“*/^T  + ‘v 

= M cosh-1  ( + N 
or 

o_  \T 
( = cosh  — ^7— 

M 

The  values  of  .If,  N,  f,  and  g can  now  be  evaluated.  At  C,  1 = 1 
and  12  = —iir/2  = N.  At  B,  t = — 1 and  12  = 0,  so  that 

7^— r=  cosh-1(—  1)  = — i cos-l(—  1)  = ±iir 

2 M 

and  M = 1/2,  the  positive  sign  being  selected.  Once  M and  N are 
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known.  / and  <i  can  bo  evaluated  by  substituting  In  (T0/T/)  and 
— i fi,  respectively,  for  fi. 


(j  = — cos  2 0 


The  equation  relating  tv  and  t can  be  obtained  as  before,  or 
directly  from  the  known  complex  function  for  the  source  and  sink 


t - plane 


t 

A 


i 


^2  ^3 

B C 


Fig.  5.  Planks  for  Schwarz-Ohrjstoffkt.  Transform atton. 


in  the  /-plane.  In  the  first  instance  the  Schwarz-Christoffel  trans- 
formation leads  to  the  form 


w=  M'J 


dt 


(<-/)  (t—g) 


■ + A7' 


The  final  form  is,  however,  obtained  without  integration  if  it  is 
noted  that  \\,b,  the  quantity  of  flow  in  the  2-planc  is  just  half  the 
strength  of  the  source  at  A and  of  the  sink  at  D in  t he  /-plane. 
Tims, 

w = — ~~  [ln(/-/)  — In  ( / — <7 ) J 


In  cases  such  as  this  one,  in  which  the  entire  boundary  in  the 
hodograph  plane  coincides  with  the  boundary  of  a circular  sector, 
a direct  transformation  from  the  £-plane  to  the  tt-plane  is  possible 
[lb].  For  the  elementary  example  herein,  this  relationship  can  be 
obtained  by  the  simple  expedient  of  constructing  the  image  pattern 
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of  the  source  and  sink  in  T ho  x-  e>"i  (/-axes  and  in  the  eireular 
boundary  ns  shown  in  Fiir.  (>.  The  vesuliimr  expression  for  w(Z) 
with  a term  for  each  of  the  four  sources  and  four  sinks  will  lead 
to  a relationship  equivalent  to  that  obtained  by  the  method  out- 
lined first.  Alternatively,  the  c-  and  /'-planes  are  simply  connected 
by  means  of  the  relationship 


Which  planes  and  which  methods  are  appropriate  ir.  a given  case 
depends  upon  the  type  of  problem  and  the  nature  of  the  calcu- 
lation to  bo  made. 

Once  ic(! } and  either  fit/)  or  £(/ ) have  been  defined,  the  desired 
function  u'(z)  is  implicitly  defined,  both  9 and  s being  expressible 
in  terms  of  du'/dz.  If  the  resulting  ordinary  differential  equations 
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can  he  integrated,  three  kinds  of  results  can  he  obtained:  (1)  para- 
metric equations  relating  buik  characteristics,  i.e..  a/b,  1\.,  and  ft, 
in  the  foregoing  example,  (:2>  the  profile  of  tin*  jet,  and  (d)  the 
variation  <>f  velocity  and  pressure  along  the  solid  boundaries.  Such 
calculations  are  presented  in  some  detail  in  each  of  the  following 
papers,  fn  addition,  a still  more  complicated  result,  which  usually 
requires  numerical  integration,  is  available  in  that  <f>  or  ^ can 
be  defined  at  any  point  in  the  region  of  flow  in  the  2-plane. 
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Arri.Kwmi.iT'v  ok  Method 

From  the  mathematical  viewpoint,  the  direct  calculation  of 
free-.si i eamlinc  flows  is  restricted  to  two-dimensional  irrotational 
flows  of  inoompressihle  fluids  which  are  free  from  gravitational 
effects  and  for  which  the  solid  boundaries  arc  straight.  As  flows  of 
real  fluids  seldom  satisfy  fui.y  these  stringent  restrictions,  the 
practical  value  of  the  results  is  open  to  question.  The  significance 
of  the  various  restrictions  can  be  shown  only  by  comparing  the 
calculated  results  with  those  obtained  from  observations  of  real 
flows. 

The  importance  of  viscous  effects,  or  of  relatively  small  depar- 
tures from  irrotational  flow,  has  already  been  studied  in  detail  for 
various  kinds  of  flow.  It  is  well  known  that  for  flows  in  which 
shear  is  not  a primary  factor,  such  as  rapidly  accelerated  flows,  the 
differences  between  comparable  ideal  and  real  fluid  flows  are  small. 

Tn  many  problems  the  importance  of  gravitational  effects  is 
already  known  to  he  definitely  secondary.  For  example,  the  coeffi- 
cient of  contraction  for  an  orifice  in  a vertical  wall  does  not  differ 
.significantly  from  that  for  one  in  a horizontal  wall,  even  though 
the  attraction  of  the  earth  affects  the  trajectories  for  the  two  eases 
in  quite  different  ways.  On  the  contrary,  for  flows  in  which  a free 
surface  has  a significant  bearing  on  the  flow  pattern,  as  in  flow 
over  a weir,  the  conventional  methods  are  no  longer  adequate. 

The  presence  of  curved  solid  boundaries  leads  to  calculations 
which  are  both  extremely  complex  in  detail  and  indirect  in  appli- 
cation. The  c-plane  is  no  longer  directly  useful,  and  the  Q-plane 
contains  lines  which  are  either  slanting  or  curved.  The  transfor- 
mation between  O and  thus  becomes  more  complicated,  even  if 
definable,  and  the  calculation  can  usually  proceed  only  if  the  form 
of  the  boundary  in  the  O-plane  is  assumed  first  and  the  correspond- 
ing physical  boundary  in  the  i-plane  is  derived  therefrom.  As  the 
latter  boundary  may  not  he  useful  or  even  significant,  a tedious 
t rial-and-error  calculation  will  he  required  for  the  solution  of  any 
given  curved  boundary. 

The  applicability  to  axisymmetrie  flows  of  the  results  computed 
for  plane  flows  can  be  assessed  only  by  the  making  of  suitable 
comparisons  as  in  references  [5]  and  [i2j  and  in  the  following 
papers.  Comparisons  have  been  made  both  between  results  of  theory 
and  experiment,  on  the  one  baud,  and  between  two  dimensional 
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au<l  axisynunetric  flows,  on  t lie  otlu'r.  Already  at  liand  is  convine- 
in«r  evidence  of  a surprisingly  close  correspondence  between  these 
two  classes  ot  flows,  at  least  insofar  as  certain  bulk  characteristics 
are  concerned.  Keen  in  cases  for  which  differences  have  been  noted, 
the  extensive  results  which  are  available  for  two-dimensional  flows 
are  useful  in  indicating  trends  to  be  expected  for  comparable 
nxisymmetric  and  three-dimensional  flows. 

Among  the  various  types  of  two-dimensional  flow  patterns 
which  can  he  analyzed  as  free-streamline  flows  are  free  efflux,  de- 
flection of  jets,  and  various  combinations  of  the  two.  Efflux  pat- 
terns have  been  calculated  for  a variety  of  boundary  geometries 
including:  those  for  infinitely  large  reservoirs  and  for  slots  at  the 
ends  of  conduits  with  various  boundary  configurations.  For  jots 
impinging  on  flat  plates,  the  characteristics  of  flow  have  been 
determined  for  finite  and  infinite  extents  of  fluid,  for  free  and 
partially  hounded  jets,  and  for  deflectors  of  various  proportions 
(each  being  composed  of  necessarily  straight-line  segments)  In 
another  type  of  calculation  the  free  streamline  has  been  used  as 
an  indication  of  the  configuration  of  cavities  in  fluid  flows. 

Additional  calculations  of  these  several  kinds  of  flows,  pre- 
sented in  the  following  papers,  include  characteristic  results  which 
are  useful  for  the  interpretation  of  real  fluid  flows.  The  calculations 
provide  information  for  specific  designs;  comparisons  between 
results  of  theory  and  of  observation  serve  to  indicate  the  appli- 
cability of  the  method. 


TRANSITION  CTRVES  OF  CONSTANT  PR EXSFRE 

l.  STREAMLINED  STRUTS 
by 

R.  Gkkhkk*  and  John  S.  McNown 


Transition  curves  along  which  the  pressure  is  an  arbitrarily 
assigned  constant  can  he  obtained  for  two-dimensional  flows  by 
means  of  an  adaptation  of  the  Ilelmholtz-Kirchhoff  theory  for  free 
streamlines.  These  curves  can  be  used  in  the  determination  of  the 
profile  of  the  upstream  portion  of  streamlined  struts  or  piers  so 
as  to  reduce  negative  pressures  and  prevent  marked  adverse  pres- 
sure gradients.  Consequently,  the  likelihood  of  occurrence  of  either 
cavitation  or  separation  can  be  significantly  diminished.  Although 
the  results  of  such  calculations  are  immediately  applicable  only  to 
two-dimensional  flow  of  incompressible  fluids,  the  scope  of  the  ap- 
plication is  less  restricted  than  it  would  appear.  The  hulk  char- 
acteristics of  certain  free-streamline  flows  have  been  found  to  be 
essentially  the  same  for  two-dimensional  and  axisymmetrie  patterns. 
Also,  the  criteria  for  reduction  of  shock-waves  in  compressible 
fluids  have  been  found  to  parallel  closely  those  for  the  reduction 
of  cavitation,  making  possib’e  the  use  of  such  results  in  compressible 
flows.  The  calculation  described  herein  was  suggested  by  Professor 
Kravtchenko  of  the  Faculte  des  Sciences  at  the  University  of 
Grenoble. 

A series  of  symmetrical  proflies  which  are  suitable  for  transi- 
tion curves,  and  which  are  at  the  same  time  readily  calculable,  can 
he  obtained  for  the  flow  pattern  shown  in  Fig.  la.  A flat  plate  CC' 
normal  to  the  direction  of  the  uniform  flow  at  infinity  is  connected 
to  linear  boundaries  such  as  DE  paralleling  the  direction  of  flow 
by  curves  such  as  CD  which  are  not  defined  geometrically  but 
along  which  the  vci  >city  is  equal  to  an  arbitrarily  selected  constant. 
This  constant  velocity  V,  is  greater  than  the  veloe-” 


•ity  at  infinity  V 


0> 

the  problem  reducing  to  the  classical  problem  of  the  infinite  cavity 
^ 1-1  j if  \\  is  equal  to  V„.  Considering  only  the  upper  half  of  the 
symmetrical  flow  pattern,  one  finds  that  the  direction  of  flow  is 
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a known  constant  along  tin*  segments  AB.  BC,  and  DE,  and  that 
along  the  remaining  segment  Cl)  ihe  velocity  magnitude  is  con- 
stant. Thus,  the  necessary  requirements  for  the  application  of  the 
free-streamline  theory  are  fnitiiicd. 

From  the  mathematical  point  of  view,  the  problem  reduces  to 
the  determination  of  the  complex  function  relating  the  velocity 
potential  and  the  stream  function  to  the  desired  flew  pattern  by 
means  of  an  intermediate  stage  which  is  a logarithmic  hodograph. 
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If  the  complex  potential  is  designated  as  \v  = <f> i 6,  the  complex 
velocity  can  be  represented  by  the  expression: 


d w 
~dz 


u — i v = F0e_!i 


U) 


in  which  n = — )n(fy/ro)  -f-  i 0 is  defined  in  terms  of  the  direction 
and  magnitude  of  the  vector  velocity.  The  inclination  of  the  veloc- 
ity to  the  ./--axis  is  0,  and  the  -clocity  magnitude  is  q.  The  bound- 
ary streamline  is  a semi-infinite  rectangle  in  the  Q-planc,  Fig.  lb, 
and,  of  course,  a straight  line  in  the  ic-plane,  Fig.  1c. 

The  transforming  function  between  tv  and  fi  can  be  obtained 
directly  by  means  of  the  Schwa »/.-( ‘hristoffol  theorem.  Afler  inte- 
gration  and  evaluation  of  constants, 


o = 


i . 
- cos-' 


m ( w+K) 
w 


cosh-1 


m 


(2) 


in  which  m is  defined  by  the  relationships  (see  Fig.  lc), 
K m K m 


b = 


w - 1 


a = 


///  + 1 


( m > 1 ) 


and  K is  a constant  having  the  same  dimensions  as  tv.  Along  the 
are  CD, 


so  that 


In 


fo 


= — cosir 


m 


III 


(3) 


Completion  of  the  analysis  consists  of  the  algebraic  definition 
of  the  transition  curves  as  a function  of  the  assigned  velocity  ratios, 
in  other  words,  x/c  and  iy/c  (Fig.  la)  are  to  be  expressed  as  func- 
tions of  F,/F„,  or  m.  1 localise  q — r = — ctfr/cij  along  BC,  the  ref- 
erence length  c can  be  defined  in  integral  form, 


C = f<hj  = - f ^ 

j n Jo  V 

From  Eq.  (2)  it  can  be  shown  that 
r = Fj  exp  j”—  3 cosh 


(4) 


m ( <t>  + K ) l 


J 


whereby  d<p  can  be  expressed  in  terms  of  v.  If,  in  addition,  the  sub- 
stitution t = Fj/r  is  made, 
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Along  the  are  CD.  t he  relationships  between  «/>.  0,  and  the  co- 
ordinates x and  //  are  well  known: 

I’,  da*  = — eos  6 <l<b 
V,  du  = —sin  0 d6 

Also,  from  Kq.  ( 2> . 


cos  20  = — • — - 

<*> 


/ m/ 

-T  -,  < <P  < 


V m -f  1 


Ml  — 1 / 


After  substitution  and  integration,  parametric  equations  for  the 
curve  are  obtained  in  the  following  form: 


x = 


niA 

2?; 


COS0 


Ml+1 


L 2 


COS-0 
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y=  c + 


miA 
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mi-1  m-1 

— 2~  + sm-0 


]tan~l  \l — --r  — tan*1  ( \f-  2 - - sin0 
Vm-1  \ mi  — 1 V 'Mi  — 1 


The  values  for  the  coordinates  of  point  7i  (X  and  Y in  Fig.  la) 
are  obtained  from  Eq.  (6)  if  the  value  of  zero  is  substituted  for  0: 


X = 


mi  A 

2Yi 


7„h-'U+i  i 


Y=  c + 


Ml  A 
2 F, 


I 2 


Ml  — 1 


tair 


Ml  — 1 
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Finally,  the  parameter  m,  which  has  been  defined  in  terms  of 
Vj/l'o,  can  be  expressed  in  terms  of  th  dimensionless  pressure 


change  a, 


a = 


Pn-  Pi 
pV o2/2 


(8) 
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which  is  equivalent  to  tiie  cavitation  number  il’  p1  equals  t lie  vapor 
pressure.  The  relationship  between  m and  c r is  as  iVno».. . 

" - sO+'+rb]  ,9) 

Results  of  calculations  based  on  Eqs.  ((>')  and  (7)  are  shown  in 
Fig.  2.  Tlie  solid  lines  indicate  the  shape  of  the  free  streamlines 
for  (j-values  of  0.4  and  0.7.  The  broken  line  is  the  locus  of  the 
downstream  end  of  the  curved  portion,  that  is,  it  indicates  the 
relationship  between  Y/c  and  X/c.  On  the  latter  curve  representa- 
tive values  of  o are  shown,  those  for  0.4  and  0.7  corresponding,  cf 
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Fig.  2.  Cmak.wti  ristic  Transition  Cfkvi.s  with  I.oces  of 

Tkrminat.  Points. 


course,  -with  the  downstream  ends  of  the  transition  curves  calcu- 
lated. It  is  known  that  both  the  length  and  the  width  of  the  cavity 
approach  infinity  as  a approaches  zero,  and  it  can  be  shown  from 
Eq.  (7  ) that.  A'  approaches  Y and  that  both  approach  zero  as  <? 
approaches  infinity.  From  the  practical  viewpoint,  however,  ihe 
occurrence  of  separation  prevents  the  attainment  of  flows  com- 
parable to  <r-valiies  greater  than  approximately  0.5,  a point  dis- 
cussed by  Birkhoff  1 13]. 

Although  marked  differences  exist  between  results  of  compu- 
tations based  on  two-dimensional  flows  and  observations  of  other- 
wise comparable  three-dimensional  flows,  remarkably  close  com- 
parisons have  been  found  between  certain  bulk  characteristics 
[0,12.15].  Whereas  the  profiles  of  free  streamlines  are  markedly 
different,  the  ratios  of  comparable  areas  have  been  found  to  be  the 
same  for  the  two  kinds  of  flow,  an  observation  which  led  the  writers 
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Fie.  3.  Comiwrimix  of  Ratios  Cai x ti.atku  for  Tv. o-Dimknsionao  Flow 
with  ntosi:  Mkasi  ki  l)  FOR  Axisymmktkh:  Fi.OWS. 


to  make  a similar  comparison  of  the  IVc-valucs  with  the  ratios  of 
t he  area  of  a plate  or  a disk  to  the  measured  cross  section  of  cavi- 
tation pockets  they  cause  to  form  in  axisymmetric  flows.  The  com- 
parison between  the  results  obtained  from  Eq.  (7)  and  the  obser- 
vations made  in  the  water  tunnels  of  the  Iowa  Institute  of  Hy- 
draulic Research  [10]  and  of  the  David  Taylor  Model  Basin  [17] 
is  shown  in  Fig.  *».  Surprisingly  good  correspondence  is  obtained 
for  small  values  of  <r;  the  effect  of  separation  is  evident  in  the 
discrepancies  occurring  for  a greater  than  0.6.  The  agreement 
obtained  is  better  than  that  which  Risen  berg  [IS]  found  in  a sim- 
ilar comparison  with  the  theory  of  Riabonchinsky  and  that  of  Oil- 
barg  and  Rock,  but  the  reason  for  this  improvement  is  not  evident. 

In  conclusion,  the  writers  believe  that  the  calculation  of  the 
irrotational  flow  patterns  herein  presented  is  directly  applicable 
in  certain  design  problems,  that  it  can  be  extended  or  modified  to 
describe  other  related  flows,  and  that  certain  applications  are  pos- 
sible for  axisymmetric  flows. 


TRANSITION  (TKVKS  OF  CONSTANT  PEESSFRE 

II.  INLETS 

by 

I).  W.  Arpr.i,  and  E.  M.  Lvuksf.n* 


Patterns  ot‘  flow  through  well  stream  line*.  .1  CO  1 1 tractions  have 
been  found  to  be  detined  at  least  moderately  well  by  potential-flow 
theory.  The  particular  type  of  contraction  transition  to  be  con- 
sidered herein  is  characterized  by  a constant  pressure  (and  there- 
fore a constant  magnitude  of  velocity)  along  the  curved  portion 
of  the  boundary.  This  typo  of  transition  for  two-dimensional  flows 
is  amenable  to  free-streamline  analysis  provided  that  the  boun- 
daries upstream  and  downstream  from  the  constant-pressure  sec- 
tion are  made  up  of  straight  lines,  along  which  the  direction  of 
the  flow  is  constant.  The  form  of  the  curved  boundary  thus 
should  lie  calculable  by  the  free-streamline  techniques  outlined  in 
the-  first  paper  of  this  bulletin  and  t lie  problem,  is  an  extension  of 
that  in  the  preceding  paper. 

Several  advantages  are  inherent  in  this  constant-pressure  transi- 
tion. The  minimum  pressure  in  the  transition,  that  along  the  free 
streamline,  can  be  specified;  thus,  the  likelihood  of  cavitation  can 
be  directly  controlled.  Since  there  is  no  adverse  pressure  gradient 
within  the  curved  portion  of  the  transition,  the  probability  of  a 
separation  of  the  flow  from  the  boundary  is  also  reduced.  Finally, 
tin*  occurrence  of  the  maximum  velocity  all  along  the  free  stream- 
line should  result  in  a minimum  size  for  the  transition. 

The  location  of  the  boundaries  of  a general  two-dimensional 
contraction  arc  shown  in  Fig.  la.  For  an  infinite  value  of  the  up- 
stream width,  the  figure  represents  one  half  of  a symmetrical  inlet 
from  an  infinite  reservoir.  For  upstream  widths  greater  than  the 
downstream,  the  figure  represents  various  degrees  of  contraction 
of  two-dimensional  flow.  If  the  centerline  and  the  upstream  boun- 
dary are  considered  to  he  repeating  lines  of  symmetry,  the  figure 
represents  either  the  upstream  portion  of  a series  of  stmts,  or, 
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ly.  a series  of  equidistant  inlets  from  an  infinite  reservoir. 
arSsw:  iv  he  derived  for  the  form  of  the  constant-pressure 
■j.  can  he  applied  in  >3»t«ie'-  rt  flow  past  any  of  these 


boundary  forms.. 
M-athehath  ail  Analysis 


The  four  j •lanes  described  in  the  introductory  paper  and  used 
in  the  following  analysis  of  two-dimensional  How  contractions  are 
shown  in  Fig.  1.  In  the  flew  plane  ■ the  r -plane).  the  inference 
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velocities:  along  the  boundary  are  shown.  The  flow  is  represented 
in  the  logarithmic  hodograph  (She  fi-p!anet  by  the  usual  relation- 
ship, 

f!  = — In  ~ + ?"  & i 1 ) 

* T> 

for  which  the  bounding  streamlines  take  the  form  of  a semi-infinite 
rectangle.  The  open  end  of  the  rectangle  at  infinity  corresponds 
to  tire  stagnation  paint  *4.  and  the  Swo  90=  vertices  correspond  to 
points  II  and  C at  each  end  of  the  free  streamline. 

By  means  of  the  Schwarz-Chridaffel  theorem  the  semi  infinite 
rectangle  in  thr  ft-plane  is  mapped  onto  one-half  of  the  auxiliary 
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i1 -plane..  Poii!t>  .1  and  /»  «<•*■  »rh:*r;  liiv  t*hu-e>3  at  zero  and  infinity 
in  the  f-plane.  Rat  her  than  to  fix  V’;  -'•Miit  in  u-  plane,  it  is 
convenient  to  relate  the  distance?  / ani  *,  by  the  paru;netri»-  eqna- 
1 ions 
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The  Seh’.varz-ChristoF*-]  transformation  between  the  £l-plane  and 
the  f -plane  is  tlieu 


o 


i;2> 


The  constants  in  the  transformation  were  evaluated  by  substitu- 


tion of  the  values  of  Q and  r for  points  A and  li  in  the  two  planes. 
An  cognation  for  the  parameter  m can  be  obtained  by  substituting; 
values  of  11  and  r for  the  point  D in  the  two  planes,  whence 


--5[(r:)+(f:n 


Thus,  if  the  velocity  alone  the  transition  V.  relative  to  the  velo- 
city in  the  downstream  channel  r>7,  is  specified,  m is  determined. 

The  infinite  rectangle  representing  t (he  potential  and  stream 
functions  in  tBte  ir-pkmv  is.  mapped  onto  the  lower  half  of  the  f-pBaiue 
by  the  transformation 

?r  = <i>  i $ — — — In 1 i k I i 3 ) 

1 T 


This  equation  can  be  derived  by  means  of  the  Schwarz-Christofifel 
theorem  or  from  the  expression  for  a source  located  at.  point  E. 
The  evaluation  of  Eq.  (2  :■  for  the  point  E gives  an  expression  for  h. 


It — TO 


in  which  >s  is  determined  by  the  ratio  of  the  free-streamline  velo- 
city 1'.  to  the  upstream  channel  velocity  T,  according  to  the  equa- 
tion 
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If  the  restriction  that  the  magnitude  of  the  variable  velocity  q 
equals  that  «.»t  » is  imposed,  parametric  differential  equations  for 
the  free-stnanJiru  can  be  obtained.  Equation  1 ? thus  becomes, 
for  the  free  streamline. 
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9 = — !n  ^ i 0 (4) 

* o 

Simultaneous  solution  of  Eqs.  (2),  (3).  and  (4)  gives  the  follow- 
ing relationship  between  the  complex  potential  ic  and  t ho  incli- 
nation of  the  free  streamline  0: 
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For  eonvenienee  the  free  streamline  has  been  taken  as  f = 0.  The 
distance  s along  the  free  streamline  from  B to  a point  at  which, 
the  inclination  is  6 can  be  obtained  by  combining  the  equation  de- 
fining the  velocity.  d<f>/cs  = —V>,  with  Eq.  (5),  and  integrating: 

V0  a . T to 
2ttV<i 


= i"  r- 

2ttV->  L ti 


+ 1 
n—vi  J 


(6) 


to  — eos  2 0 

Because  s is  the  curvilinear  distance  along  the  bounding  stream- 
line, this  expression  is  not  directly  useful  for  obtaining  the  profile 
of  the  transition.  However,  the  equation  can  he  expressed  in  rect- 
angular coordinates  through  the  differential  relationships 

dx  = ds  eos  v . dy  — <»<?  sin  0 

The  parametric  differential  equations  for  the  curved  boundary  are 
then 
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. ( to  — 1 + 2 sin20 ) (n  — i.  + 2 sin20 ) . 

Fortunately  these  equations  are  readily  integrable;  one  obtains 
the  parametric  equations 


x l V0  r * i , cos0 

o = -vr[v!».+  a/2 


COS0 


' (»+n/2ta»b->  ®> 


and 

y 

a 


sin0 

l)/2 


1 , 1 Vo  r /7- — t -i i?1 

2 + 7rF2LVU  1 ) — v/(n 

-V(^PTPt.n-.7jMT7l]  (10) 


Transition  Cirvks  of  Constant  Prkssvrf — II 


if  t ]i r>  origin  of  ./•  and  >i  :s  chosen  such  that  x = 0 if  8 = — r/2  and 
7 =-  a 2 if  0 = 0.  Equation  ( 9 \ can  he  written  in  an  alternative 
logarithmic  form. 
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For  the  case  of  the  infinite  reservoir  in  which  n = oc,  Eqs.  (9) 
and  (10)  can  be  reduced  to 
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after  evaluation  of  the  indeterminate  forms.  Rethy  in  1893  |19] 
presented  equations  similar  in  form  to  Eqs.  (9)  and  (10).  How- 
ever. his  equations  were  based  upon  parameters  which  could  not 
he  directly  evaluated  from  the  flow  characteristics. 

Upstream  and  downstream  from  the  curved  boundary  the  magni- 
tude of  the  velocity  is  variable,  hut  the  direction  of  the  velocity 
along  the  boundary  is  known.  Equations  (1),  (2),  and  (3)  can  be 
used  to  obtain  the  velocity  distribution  along  the  straight  portions 
of  the  boundary  by  steps  essentially  parallel  to  those  used  for  the 
calculation  of  the  free-streamline  profile.  For  t lie  general  case,  the 
velocity  distribution  can  be  obtained  from  the  equation 


+ K4  (tanh-'  ^ -Unh-1^)]  (13) 

downstream  from  the  transition,  and 
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upstream  from  the  transition.  In  Eqs.  (13)  and  (II), 

Kx=  [(F2/F0)2-  (F./Vjii-t  m/vny- (Ki/n)2]-1 
K2=(v0/v2y  [(Vo iv2y-{v2iv,y)-'  {(vjv2y- (Fi/f»)*]-» 
KZ=(VJ\\)  [(F2/Fi)2  — (V'i/Fo)2]-1  [(V2/Fi)2—  (F0/F2)2]-1 
Kr^iVoVJVS)  [(F,/F2)2-(K2/To)2]-1  (OVF,)2-  (Fo/F,)2]-1 


and  c and  d are  the  dimensions  of  the  transition.  The  form  of  these 
equations  for  the  special  ease  of  the  infinite  reservoir  is 


a a 7r  r 2 L q J x \K2/  L K2Z  F2J 

-f-  Ftanh-1  — — tanh-1  Fi  1 (15) 

x L g K2  J 
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upstream  ftom  the  transition.  The  pressure,  distribution  can  be 
obtained  directly  by  means  of  the  Bernoulli  theorem. 

Particular  constant -pressure  transitions  can  now  be  calculated 
for  specified  values  of  the  velocity  ratios  YJY0  and  V.,/Vly  or, 
what  is  equivalent,  for  specified  values  of  the  minimum  pressure 
along  the  transition  and  the  contraction  ratio  b/a.  In  order  to  make 
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3he  results  of  sample  calculations  wncral.  values  of  a dimension- 
less pressure  parameter  <x  have  been  assumed  and  the  correspond- 
ing pressures  calculated  from  the  definition  equation. 
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Pi  ~ P 2 


(17) 


If  p2  is  the  vapor  p~  ssure.  s is  the  familiar  cavitation  index.  Since 
the  pressure  and  velocity  an?’  related  by  the  Bernoulli  equation,  toe 
velocity  ratio  1*.  Yh  can  also  be-  obtained  from  the  equivalent  ex- 
pression 
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Transition  profiles  have  been  calculated  for  a value  of  <?  of  0.4 
and  contraction  ratios  b *i  of  ac»  4.  artel  2 The  profiles  are  shown 


o 

Pks.  2.  Tr«sniw  Profiles  ram  a Constats  t Petsscei-  Ixmx  or  <r  = 0.4. 

in  Fig.  2.  For  different  values  of  «r.  the  two  principal  dimensions 
of  the  transition,  the  lengtn  d and  the  width  c,  are  shown  in  Fig.  3. 
As  a approaches  zero,  i.e.,  as  1\ — *1',,.  the  length  of  the  transition 
approaches  infinity. 
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Comparison  UiinvEEX  Results  of  Thjmry  and  Experiment 

Iti  1917.  vom  Misers  !4]  published  an  expression  lor  tire  con- 
traction coefficient  for  free  efflux  corresponding  to  V in  the 
above  analysis.  Equation  ( 10  for  m = 1 anti  9 = — sr  2.  iff  ex- 
pressed in  terms  of  a contraction  coefficient,  is  found  to  be  iden- 
tical to  von  Msses’  equation: 
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Although  von  Mks  did  not  compute  nhe  jet  profile,  he  used  the 
free-streamline  theory  to  obtain  the  pressure  distribution  on  the 
channel!  walk  With  knowledge  of  the  {pressure  distribution,  lie 
obtained  the  contraction  coefficient  by  means  of  the  equation  of 
moment  tint.  Since  von  Miss’  equation  is  a limit  to  the  theoretical 
expressions  obtained  here  so.  observed  agreement  between  the  re- 
sults of  ex|erimp«5s  and  his  computed  coefficients  is  a partial 
verification  of  the  present  analysis. 

Ira  order  to  test  the  performance  of  the  theoretical  transition 
forms  in  detail,  experiments  were  coudm-ted  with  air  flow  through 
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hah  nK«deb  >*J  «*r‘  th»  prifilrs  shown  in  Fi’».  2 {b, <?  = ac  and 

h ’**  = — . 2 lie  curved  bot  inflow  .iml  rv»rt  ions  of  tiss?  st  •rsnsiit  botus- 
dary  upstream  rut  \ downstream  were  made  of  Lueite  with  r<iezo- 
m**l  r:e  taps  along  1 j'j >*  centerline.  The  width  «i  ‘2  was  9.2  inches  and 
th*1  lictsrftt  w;js  also  >.2  inches.  makins  the  downstream  cross  section 
square.  Tb*  flow  of  air  was  produced  and  controlled  by  a variable* 
disclrftrgt-  blower  located  seven  feel  downstream  from  the  transition 
seer  son.  The  rate  of  llw  was  determined  by  means  of  a calibrated 
:joz£.iie  os?  tin-  diseloao ge  side  or  the  blower. 

The  r li eoretie.il  pressure  distribution  along  the  boundary  of  the 
in  let  frem  a he  tniinire  reservoir  as  shown  in  Fig.  4a  along  with 
test  results  for  downstream  conduit  velocities  of  IT  and  46  ffps, 
corresponding  to  Reynolds  numbers  of  To.000  and  200,000  based 
ot3  the  halt-width  of  9.2  inches.  An  approximately  constant  pres- 
sure was  obtained  along  the  curved  boundary  for  both  tests.  The 
greater  reduction  in  pressure  for  the  higher  velocity  indicates  a 
trend  toward  the  theoretical  value  for  ct  = 0. 1 with  increasing 
Reynolds  number.  That  viscous  effects  result  in  a higher  pressure 
is  note  wort  by.  because  the  design  thus  tends  to  be  conservative 


Fig.  I.  I’bskssl-kf.  Bi-shtriritioxs  Ovtaixkd  Along  Bocndaeiks  or 
Ponvro  Tb-ixssI'bon  Fiwses. 
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and  the  likelihood  of  cavitation  tends  to  be  even  less  than  antici- 
pated. 

The  pressure  distribution  along  the  boundary  for  the  2:1  chan- 
nel contraction  is  shown  in  Fig..  4fe  for  the  same  Reynolds  numbers 
a*v  above.  The  velocii*  distribution  in  the  upstream  channel  evi- 
dently affects  the  pressure  distribution  .so  that  the  pressure  along 
the  transition  boundary  is  net  constant  and  the  minimum  pressure 
occurs  near  the  downstream  end  of  the  transition.  Again  the  trend 
is  to  approach  the  theoretical  pressure  reduction  at  higher  Reynolds 
numbers.  However,  the  eddy  which  forms  near  the  right  angle  at 
the  beginning  of  !he  contraction  would  probably  affect  the  pres- 
sure distribution  regardless  of  the  Reynolds  number. 

Comparison  with  Other  Transition  Forms 

Figure  4a  also  shows  the  results  of  a lest  that  was  made  on 
a nr  odd  off  one  of  the  constant-pressure  transitions  In  an  electrical- 
analogy  tank.  Since  the  electrical-analogy  method  is  based  on  the 
same  assumption  of  potential  flow  as  is  implicit  in  the  free-st  ream- 
lire  analysis,  the  excellent  agreement  with  the  theoretical  pressure 
distribution  simply  verities  the  technique  and  accuracy  of  the 
electrical-analogy  method.  However,  the  check  allows  the  results 
of  tests  (as  yen  unpublished)  on  circular  and  elliptical  transitions 
to  be  used  with  confidence  in  the  comparison  of  transition  dimen- 
sions shown  in  Fig.  b.  Although  the  data  available  from  the  elec- 
trical-analogy experiments  are  limited  to  a few  shapes  and  to 
inlets  from  an  infinite  reservoir,  the  surmise  tliat  die  dimensions 
of  the  constant -pressure  transition  arc  the  smallest  possible  for  a 
given  minimum  pressure  along  the  boundary  is  upheld. 

In  other  studi  s certain  bulk  characteristics  of  axisymmctric 
flow  have  been  found  to  be  essentially  the  same  as  those  of  two- 
dimensional  free-streamline  flow.  Two  additional  electrical-analogy 
tests  were  run  to  investigate  the  theoretical  performance  off  the 
two-dimensional  free -streamline  profile  as  an  axisymmetric  transi- 
tion. In  one  test  the  calculated  profile  for  o = 0.4  was  usrd.  i.c., 
r/ra  = y/(a/2).  The  dimensionless  pressure  parameter  at  the  be- 
ginning of  the  curved  boundary  was  positive  and  decreased  to 
about  $4  of  the  theoretical  value  of  —0.4  close  to  the  end  off  the 
transit  ion.  The  higher  pressures  were  probably  due  in  part  to  the 
larger  relative  area  at  each  section.  In  the  second  test  the  profile 
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was  modified  J«  have  the  sanx*  relative  area  at  any  spftinn,  i.e..« 


y 

'~cf  \ tj  2 

A very  large  negative  pressure,  four  times  the  theoretical,  resulted 
at  the  beginning  of  the  curve,  with  a rise  to  about  one-half  the 
theoretical  in  She  last  half  of  the  transition.  The  large  negative 
pressure  at  She  beginning  of  the  transition  is  probably  due  to  the 
considerable  curvature  at  this  point.  That  both  the  area  of  flow 
and  the  boundary  curvature  are  important  in  determining  The 
pressure  along  She  boundary  is  apparent  from  these  results. 
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•roS!?  r-tu-esentiw"  a compromise  bs'tuwn  the  eqnai- 
H|uai-iiird  profiles  tested  is  m^-essary  to  giv,  son- 
in  axbyniinetrie  flow.  Since  a mathematical  toeh- 


c • pje  comparable  «*»  eonformal  tr&nsfor  mat  ions  is  »ol  available, 
r h ■ general  solution  of  the  axisymmetrie  transition  cannot  be  ol>- 
-ti.  .Specific  csse>.  however.  earn  be  treated  by  relaxation  tech- 


UKpK.'i. 


Co>- Cl . d'.S  ION'S 

Equations  have  Been  obtained  for  the  profiles  of  constant-pres- 
sure two-dimensional  transitions  and  for  the  pressure  distribution 
upstream  and  downstream  from  the  transition.  Perhaps  the  greatest 
advantage  of  this  type  of  transition  is  that  the  minim;™  pressure 
can  he  specified.  A comparison  with  other  forms  supports  the  sup- 
position that  the  const  ant  -pressure  transition  is  the  smallest  pos- 
sible for  a given  minimum  pressure  along  the  transition. 

The  limited  experimeivts  with  flowing  air  indicated  good  cor- 
respondence between  results  of  theory  and  of  real  fluid  flows  for 
high  Reynolds  numbers  ami  for  large  contraction  ratios.  For  lower 
Reynolds  numbers  and  for  smaller  contraction  ratios  the  corres- 
pondence was  not  as  good  because  of  the  effects  of  viscosity. 
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DEFLECTION  OF  JET.*> 

SV\!\IETRK'ALLV  PLACED  V-SHAPED  ORSTAIkE 

by 

T.  T.  Si ao  and  P.  G.  Hubbard* 


The  d-r  fleet  Son  of  a free  jel  In*  a solid  boundary.  which  has  Ions: 
fee 3 1 utilized  to  develop  power  from  flowing  water,  is  well  suited 
to  free-streamline  analysis  because  of  the  dominance  of  inertia  and 
pressure  intensity  in  the  establishment  of  the  flow  pattern.  The 
design  of  impulse  machinery  utilizing  this  momentum  ehange  could 
he  facilitated  greatly  if  She  idealized  geometry  of  the  system  under 
potential  flow  conditions  were  known,  because  such  conditions 
represent  asymptotic  values  which  are  approached  as  the  effects 
of  secondary  variables  ace  decreased.  With  such  information  avail- 
able. refinements  of  design  com  Id  be  based  upon  a secure  knowledge 
of  tin  f (neda mentals,  and  many  rules  of  thumb  could  be  replaced 
with  precise  quantitative  data  in  graphical  or  tabular  ferm.  Spc- 
eifieally.  if  the  total  angle  through  which  the  let  is  deflected  is 
determined  for  conditions  of  both  partial  and  complete  interception 
hy  the  itouudary.  then  the  principle  of  impulse  and  momentum  can 
he  used  to  compute  forces  or  other  dynamic  characteristics  of  the 


This  paper  and  the  two  u'hiclj  follow  are  devoted  to  a determi- 
nation of  the  angles  of  deflection  caused  by  certain  idealized  forms 
of  solid  boundaries  placed  cither  symmetrically  or  asymmetrically 
with  respect  to  the  axis  of  a two-dimensional  jet.  These  patterns  of 
flow  earn  p-snd  to  tint's,  occurring  as  a bucket  of  an  impulse 
machine  passes  through  a circular  jet.  In  this  paper,  the  free- 
streamline  method  is  used  to  find  the  angle  through  which  a two- 
dimensional  free  jet  will  be  deflected  by  a y\  iridetrical  V-shaped 
boundary  placed  on  its  axis. 

As  iT-presented  in  Fig.  la  ft  he  --plane;1.  a two-dimensional  jet 
with  velocity  1,,  and  width  2a  is  deflected  5 Si  rough  an  angle  0 by 

*Researcti  Associate  and  Research  Enviaeer.  respectively,  Iowa  Instu- 
rnif  ©£  Hydraulic-  Research.  Slat...  University  of  Iowa.  Iowa  City. 
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the  annular  boundary  of  projected  width  2b.  T7-e  sides  of  the 
b*»unda»’y  arc  inclined  at  an  angle  a relative  to  the  original  jet 
direction.  In  the  hodbgraph  plane  Fie.  l»>t.  the  bounding  slreatn- 


0 I M,  -I  0 l 

(b)  5 -plane  (c)  J'-plone 
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(d)  ? - p«cne  (e)  w -plane 

Fig.  1.  TkansvFC'Bmatios  Flasks. 


lines  become  a circular  sector  of  angle  a as  can  be  shown  from  the 
definition  relationship 


in  which  tr  is  the  complex  potential.  (Details  oi  this  concept  are 
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given  in  the  first  paper  in  this  Bulletin.)  This  sector  is  expanded 
to  a semicircle  (Fig  !<•>  by  the  transformation 


(2) 


and  finally  into  the  usual  /-plane  (Fig.  Id)  by  the  transformation 

1 “ \ (f'  + p)  (3) 

In  the  /-plane,  the  flow  pattern  is  simply  that  for  a source  with 
strength  aV j at  A and  a sink  of  equal  strength  at  7)  for  the  floiv 
occurring  in  the  upper  half  of  the  plane.  This  latter  requirement 
is  satisfied  if  the  strengths  of  the  source  and  sink  are  doubled  to 
provide  an  image  pattern  in  the  lower  half  (which  actually  repre- 
sents the  omitted  half  of  the  original  system).  Application  of  the 
usual  equations  for  a source  and  a sink  on  the  real  axis  then  results 
in  the  potential  function  (Fig.  le) 

w — Jjn(/-f-cos^~^  — ln(Z-f  1)J  (4) 

From  the  purely  analytical  viewpoint,  Eqs.  (1),  (2),  (3),  and 
(4)  represent  the  solution  to  the  problem,  because  proper  manipu- 
lation of  the  variables  will  yield  values  of  the  velocity  at  any  point 
in  the  physical  plane.  The  explicit  solution  of  these  simultaneous 
algebraic  and  differential  equations  as  they  stand  is  cpiite  involved, 
however,  and  the  complexity  of  detail  can  be  mitigated  consider- 
ably by  the  introduction  of  an  auxiliary  variable  t,  defined  by 

T = jM/n  (5) 

and  two  positive  integers  m and  n (n  < m)  such  that 
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Then,  from  Eq.  (1) 

_ I dw  _ 1 dw  dt  df'  df  , 
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and  substitution  of  the  appropriate  derivatives  from  Eqs.  (2-5) 
yields 
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The  resolving  *•;  ths-  n nation  into  partial  fraetiwis  results  its 
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Integration  of  this  equation  between  proper  limits  will  give  the 
coordinates  iti  i iti  z plane  ot  any  jiomi  corresponding  to  assigned 
values  of  a and  $ and  of  t.  The  differences  in  the  coordinates  of 
points  B and  for  example,  can  be  found  by  noting  that,  for  the 
point*  B ami  C, 

l B = 0 


or 


-#  = Q 


if  = f 


~c  = e! 


Thus,  since  B is  the  origin  in  the  ;-pl»m\  integration  of  Eq.  (6  - 
between  the  limits  of  0 and  eia  K gives:  tBie  coordinates  of  point  (*-. 
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The  imaginary  part  of  this  equation  is  the  {/-coordinate  of  point  C. 
which  is  equal  in  absolute  value  to  the  distance  denoted  as  b in  the 
s-plane: 
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Vertical  bars  in  Bq„  (7  i;  indicate  that  ahsnluto  values  are  to  be  used. 

Although  analysts  haw  already  presented  geneva!  solutions  of 
jet  interceptions,  and  Visotiti  *20*  has  completely  solved  the  par- 
ticular ease  of  the  normal  plate,  the  method  presented  foencin  is 
considered  more  direct  and  the  integrated  solution  is  completed 
for  the  general  case  in  Eq.  (7  . From  this  equation,  values  of  b/a 
corresponding  to  assumed  values  of  0 have  been  computed  for 
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various  vaiut>  of  «.  The  2 vaults  are  plotted  in  Figs.  - and  3,  and 
•n.  -vrsd  items  ♦»:'  interest  are  Lmm^liatelY  apparent.  For  example, 
a boundary  with  a projected  width  of  {trice  the  jet  tfcickue**  «ib- 
r j«*.r  <r-.i  dcih  rtioits  greater  .kata  !MF.  rvgardless  of  the  angle  at 
whirl  it  j>  plated.  Similar  conclusions  could.  also  be  drawn  for 
jst her  relative  width  ra;:**s..  < hi  the  other  lutnd,  if  the  piojtcled 
width  «*2"  the  biiiidatr  i>  four  times  the  jet  thickness,  then  the 


b/a 

Fic.  2..  KiTtrr  «*”  I>m ’.\i>.%ky  Wmac  ox  J?..t  Dimawx. 

jet  deflect  doth  will  be  nearly  complete  for  alt  angles.  This 

means.  of  that  further  widening  of  the  deflector  will  increase 

the  force  exerted  by  only  a small  amount;  "he  amount  of  possible 
increase  is  indicated  by  the  45-dcgrce  line  corresponding  to  an 
inliniteSy  wide  deflector  (Fig.  3 . 

In  order  to  compare  the  results  of  this  calculation  with  those 
observed  for  real  axasy  inmetric  flows.  cxj»ertmc«ts  were  conducted 
it;  which  a high-velocity  circular  jet  of  water  was  deflected  by  cir- 
cular disks  of  various  diameters  (corresponding  to  a — *■/ 2). 
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Larger  values  of  a would  be  used  m an  impulse  machine,  of  course, 
hut  the  normal  disk  is  readily  made  and  causes  deflections  whir  Si 
can  lie  measured  with  an  accuracy  sufficient  for  comparison  with 
computed  values.  The  deflection  jti  caused  by  each  of  three  disks 
was  measured  from  photographs  taken  with  a camera  mounted 
vertically  above  the  disk  to  minimize  errors  due  to  the  effect  of 
gravity  upon  the  stream 

Be  cause  the  phenomenon  under  study  is  primarily  one  of 
momentum  change,  and  beca  use  the  area  rather  than  the  width 
enters  linearly  into  momentum  espial  ions,  it  seems  plausible  that 
the  area  ratio  should  have  more  significance  than  the  diameter 
ratio  in  comparing  the  a visym metric  experimental  system  with  the 
two-dimensional  analytic  system.  In  the  latter,  of  course,  the  areal 
and  linear  ratios  are  identical,  in  aceordain-e  with  this  concept, 
measured  values  of  the  deflect  ion  d have  been  included  in  Fig  2, 
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the-  ratio  of  disk  awa  to  jet  cr>^s-sectionaL  area  having  been  used 
as  nBac-  abscissa.  From  these  whits..  it  can  :k-  concluded  that  the 
deflector  is  less  effective  than  computed  values  would  indi-nto.  and 
that  then.-  is  an  increasing  discrepancy  between  computed  and 
measured  values  as  the  ana  of  the  deflector  increases.  Comparative 
values  of  the  momentum  chance  are  06.  90.  and  S3  percent  of  the 
predictions  for  area  ratios  *>f  0.5.  1.  and  2.  respectively. 

The  correspondence  between  the  measured  values  and  the  com- 
pufed  ones  is  illustrative  of  the  degree  of  correlation  between  the 
idealized  two-dimensional  flows  and  comparable  real  flows  of  prac- 
tical significance.  As  already  indicated  in  the  first  paper  of  ibis 
series,  the  reasons  for  a lack  of  complete  agreement  can  bo  deter- 
mined by  finding:  the  ways  in  which  the  experimental  system  fails 
to  meet  the  assumptions  made  in  the  mathematical  analysis,  In 
addition  to  the  probable  modifications  of  the  flow  pattern  intro- 
duced by  using  an  axisym metric  system  rather  than  one  of  two 
dimensions  only,  it  can  be  expected  that  viscous,  gravitational,  and 
possibly  capillary  effects  may  modify  the  flow  to  a measurable 
degree.  Although  the  magnitude  of  the  effects  of  these  ignored 
parameters  is  difficult  to  determine,  it  can  be  shown  that  their 
influence  in  each  case  is  in  tine  projwr  direction  to  help  account  for 
the  discrepancy  actually  found.  For  example,  curvature  of  the 
pa tli  lines  of  the  deflected  sheet  due  to  gravitational  acceleration 
will  decrease  the  apparent  angle  of  deflection,  and  extraction  of 
energy  from  the  system  by  viscous  friction  at  the  deflector  and  at 
the  large  water-air  interface  will  decelerate  the  flow  and  conse- 
quently decrease-  the  angb  of  deflection.  Capillary  ciitraMf-ition  at 
the  solid-liquid-gas  juncture  where  separation  occurs  should  be 
small,  but  nonetheless  will  act  to  bend  the  flow  in  the  direction 
measured,  and  surface  tension  in  the  deflected  conical  sheet  will 
fend  to  resist  the  necessary  increase  in  total  area  with  distance  as 
required  by  continuity  principles  applied  to  the  axisym  met  trie  case. 
This  latter  effect  in  particular  would  increase  rapidly  as  the  de- 
flection angle  becomes  greater,  and  may  account  for  the  larger 
discrepancy  found  at  larger  angles. 

In  conclusion,  the  in  format  ion  made  available  ami  the  close 
correlation  between  the  results  of  free-streamline  analysis  and 
those  of  experiment  are  enough  to  warrant  tlse  recommendation 
of  this  method  as  a valuable  supplement  to  existing  methods  em- 
ployed in  engineering  design.  Care  mud.  of  course,  be  exercised 
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in  interpreting  the*  results  of  eompntat  ions  of  a two-dimensional 
flow  pattern  for  a real  system  which  is  r.xisymmetric.  as  illustrated 
by  the  use  of  areal  ratios  rather  than  linear  ratios  in  comparing 
the  present  results.  For  a system  in  which  the  Sow  parserii  is 
determined  primarily  by  inertia  and  pressure  intensity,  a knowl- 
edge of  the  true  potential  conditions  is  a valuable  asset  in  assessing 
the  effects  of  secondary  variables  and  la  cnee  in  judging  what  im- 
provements are  possible  through  iwluction  of  .secondary  effects, 
in  relation  to  the  problem  of  Pel  ton -wheel  design.  it  has  been  shown 
that  the  angle  through  which  the  jet  is  deflected  can  toe  substan- 
tially less  than  the  angle  on*  the  bucket  lip.  approaching  the  latter 
value  only  if  the  projected  area  of  the  bucket  is  at  least  four  times 
the  jet  area.  This  means  that  the  bucket  should  3*e  quite  deep  for 
practical  angles  of  deflect  ion.  and  that  the  bucket  spacing  may  be 
affected.  A combination  of  the  results  presented  herein  and  those 
of  the  following-  two  papers  may  prove  valuable  for  further  study 
of  the  Pelton-wheel  design. 
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In  the  preceding  analysis  the  problem  of  a jet  striking  a wedge- 
shaped  obstacle  was  solved  by  the  method  of  conformal  mapping. 
In  that  analysis  the  jet  in  consideration  was  free  in  the  sense  that 
it  was  not  confined  by  any  rigid  boundary  before  it  struck  the 
object.  However,  in  some  applications  of  the  preceding  calculation, 
the  approaching  jet  os  not  everywhere  free  but  is  rather  confined 
in  a pope  or  channel.  With  the  introduction  of  a straight  boundary 
Ilf'.  Fig.  4.  the  jet  is  no  longer  entirely  free  and  the  solution  must 
be  modified  accordingly. 

in  Fig.  4.  a confined  jet  of  half  width  h with  an  initial  velocity 
V , ptmeds  towards  the  right,  meets  the  harrier  AK  having  an  .ingle 
of  inclination  « to  the  dir*  et ton  •>.»!  the  approaching  jet.  and  deflects 
through  an  angle  8.  The  par  tern  obtained  by  revolving  the  figure 
about  the  axis  AH  eorresjionds  to  that  of  an  axially  symmetric  flow 
through  a cone  valve  21  . The  discharge  coefficient  is  a function 
of  the  angle  a.  the  opening  s.  and  the  dimension  S. 

Application  of  Eqs.  * 1>  and  '2  transforms  the  flow  in  the 
'-plane  into  that  within  a half  circle  in  the  f-plane,  with  a source 
of  strength  2d!’,  at  jx>5uit  It  and  a sink  of  equal  strength  at  point.  D. 
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in  whieh  c = d h.  By  virtue  of  Eq>.  f 1;.,  {2],  and  (8), 

I f 5 v 1 f 1 d sr  d f v 
~ = f ^ d nr  = - - | — - — 5-  *$c 
* } J 1 V n j % <t  i fJ  \ 

(I — 1 

+1 ’ L. ,lr...-=rff 

iT  T **  — f ^ os f ~ * 


If  she-  Mtbsflilutions 


42 


Fes  f -5nr../\'o;M ve;  Asai.vsi.  > 


n 

Ct  — X 


an*  asa’m  niade.  m and  » being  any  integers  provided  lliat  .?»  > n, 
Hu  i;lfj  reduces  7© 


I 


T 


— c 


4- 


1 

— ^ £ — S T 


£ y*L  f i _ ]e  — -ra C 

- > c i1  .. rf-l  *4  !:.?«•  **•  ? 

»•  ♦.  • * *L  A tL  *1- 

y = i? 

f~  » ,T  -iCirar—  5-  jf 

_jl | ; •.  i _ 

■'  <v  ' j-.  > • I * ffl  — S* ’ m | ** 


(10) 


The  integration  of  Eq.  •>  10t  carriers  out  from  - = 0 to  r = 1 or  Ho 
r = e gives  vat ucs  of  z.  or  z itspeethvlv.  The  real  part  of  *_ 
is  equal  ito  b wt « —5.  ard  the  imaginary  par*  of  2 as  equal  to 
8-t-b.  The  final  results  are  two  equations  expressing  s and  $, 
respectively.  as  f millions  of  r,  a.  and  $. 

For  the  special  ease  of  a = x/4.  m is  4 and  a is  3 ; the  final 
results  a ire 
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and 


^4  = 2 cos  3 In  tan  ^ + 2 sin  3 In  tan  - ( ^ — 3^ 
c 22  \2  / 

-Far (cos 3 + sin 5 ) +2(  - 4-  c^tanh-3 c Hr 2{ - ; — c ^fan-3 c (12) 

the  values  of  the  arctangent  being  taken  between  0 and  x.  2.  As 
s — *■  oc , c — * 1 . and  from  Eqs.  i ll  i and  (12),  one  obtains 

— ‘ cos  3 — sin  3)  In  tan  4{t  + 

4-  (cos  3 4*  sin  3)  in  tan  ^ (^  — 3y4-  rcosd— 2 Ir-tan^  (13) 
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which  eorrtijKKwls  with  what  has  been  (kriveii  before  without  the 

boundary  Ui\  Fer  ai#»l  ’it  r va'.ivuk-  case.  .»•  — * 0 ami  < — » 0.  all  dimen- 
sions must  hr  expressed  ir  ratios  i*j  s.  Equations  <11  and  02) 
reduce  te- 
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The  above  analysis  constitutes  the  preliminary  investigation 
leading  to  the  discussion  by  one  of  the  writers,  of  the  paper  “Char- 
acteristics of  Fixed  D^ersioii  Cone  Valves'*  ] 21 1 . In  that  dis- 
cussion the  evaluation  of  the  discharge  mfficiciit  for  cone  valves 
by  an  analysis  of  the  corresponding'  plane  flow  and  a comparison 
of  the  calculated  insults  with  the  test  ones  are  given  in  detail;  in 
the  following,  only  its  essential  points  are  given. 

The  discharge  coefficient  is  obtained  by  dividing  the  discharge 
by  the  product  of  the  flow  area  and  \'%//.  Put  in  equation  form, 

C = - ® — (16) 

A\  2 gfl 

in  which  .1  is  the  area  of  efflux  and  II  is  the  total  head.  For  the  two- 
dimensional  flow  the  corresponding  discharge  coefficient  is  d 

A logical  approach  to  the  adaptation  of  these  results  for  plane 
flow  to  the  determination  of  the  discharge  coefficient  for  the  eor- 
respomiing  three-dimensional  flow  is  to  assume  that  corresponding 
ratios  between  the  initial  and  the  final  area  of  the  jet  are  the  same 
in  each  case  < ns  they  were  found  to  be  for  the  orifice  , TV  >,  That  is. 

,,  d tin-  ultimate  value  of  iri2;rr.< 

c = = — — — a it 

•*  s ( 2tBB  ) 

in  which  sc  is  the  thickness  of  the  three-climensional  jvt  at  a certain 
point  ami  r the  distance  from  that  point  to  the  g-axis. 

The  values  of  f and  & have  ?:<ecia  plotted  in  Fig..  3 against  I). 
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(which  is*  equal  to  „*  "2&  i.  A moon  curve  has  Wii  drawn  for  the 
test  results  obtained  by  Messrs.  Elder  au«l  Dougherty.  The  cor- 
resj.torodenee-  Between  a lie  text  results  and  those  of  the  analysis  ns 
fairly  good  for  values  of  v />  no  to  about  0.3. 

The  present!  |uv»Wetn  has  been  t reate*.!  by  G milfoil]  [22.  j*j». 
26-29[  who  gave  the  essential  steps  of  the  solution.  However,  the 
solut ion  pax-sentled  in  the  foregoing  is  oameh  simyder  and  has  Sxen 
carried  to  the  u tunic  rival  stage. 


DEFLECTION  OF  JETS 


JL  S V MM  E T I\T  C A LL  Y PLACED  (J-SHAPED  OBSTACLE 

by 

Turgut  Sarpkava* 

As  an  alternative  to  the  approximation  of  a bucket  on  a Pel  ton 
wheel  in  the  preceding  paper,  a channel-shaped  dellecior  can  be 
considered.  As  pointed  out  by  Greenhill  [21],  a shape  for  which 
the  irrotfltional  flow  pattern  can  be  computed  is  a symmetrical 
plate  with  two  right-angled  bends  which  would  be  tangent  to  the 
usual  bucket  as  shown  in  Fig.  1.  This  approximation  in  conjunc- 
tion with  that  presented  in  the  preceding  paper  should  provide  an 
insight  into  the  relative  importance  of  the  bucket  geometry  in  a 
determination  of  the  angle  through  which  the  jet  is  deflected. 

As  indicated  in  Fig.  2a,  a jet  of  width  2d  impinges  on  a U-shap- 
ed boundary  with  web  height  2 a and  flange  width  b,  and  each  half 
of  the  jet  is  deflected  ih rough  an  miglo  r — ft  From  the  Bernoulli 


* '{nation  the  velocity  of  approach  F,  tbe  incident  jet  is  seen  to 
be  the  same  as  the  asymptotic  velocity  of  the  deflected  jets,  and 
from  t lie  equation  of  continuity  the  thickness  of  each  of  the  de- 

♦Research  Assistant,  Iowa  Institute  of  Hydraulic  Research,  State 
University  of  Iowa.  Iowa  City. 
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fleeted  jets  is  seen  to  be  precisely  d The  total  force  exerted  by 
the  fluid  on  the  .solid  boundary  can  he  computed  from  the  mo- 
mentum principle  once  the  ancrle  of  inclination  of  the  deflected 
,jeN.  which  depends  only  on  n/d  and  h d,  is  known.  Thus,  defini- 
tion of  the  relationship 

3 = f (a I d , b ! d)  (1, 
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(o)  z-plane 


(b)  dl-plane 
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( c ) w - plane 


(d)  t-plane 


Fig.  2.  Transformation  Flanks. 
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in  the  tn;i in  objective  of  the  present  nvestigation.  For  convenience, 
the  point  of  maximum  velocity  a to  r the  solid  Imindary  between 
E and  -1  is  designated  by  the  Letter  F. 

If.  as  in  the  other  papers  eontaire.-i  in  this  Bulletin.  : and  sc 
are  used  for  the  complex  » .triable  and  the  complex  potential,  then 
the  variable  s defined  as  dsv  dz  is  given  by 

f = —u  -Eii'  = —qe~ li  (2) 

in  which  w and  r aw  the  velocity  components  in  the  x-  and  ^-direc- 
tions. q is  the  magnitude  of  the  velocity  vector,  and  0 is  its  angle 
«f  inclination. 

If  the  variable 

Q = In  \r)  = — In  0 >3* 


is  introduced,  the  flow  region  in  the  physical  plane  { Fig.  2a)  can 
be  mapped  onto  a corresponding:  region  in  the  Q-plane,  as  indicated 
in  Fig.  2b.  This  latter  region  is  in  turn  mapped  onto  the  upper  half 
of  the  auxiliary  t'-plane. 

As  indicated  in  Fig.  2d.  the  values  — 1,  0,  and  1 are  assigned 
to  the  points  D , F.  and  H in  the  e-plane.  The  values  of  f for  Er  A 
and  ('  can  he  denoted  by  —k.  k,  and  /.  respectively,  t loose  for  E 
and  .1  being  of  equal  magnitude  and  opp«3site  signs,  due  to  sym- 
metry. The*  desired  mapping  between  U and  i is  Inen  achieved  by 
the  use  of  the  Schwarz-Chiistolnfe!  transformation 
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in  which  A is  found  to  be  \ 1 --  k*  by  substituting  into  the  equa- 
tion  the  coordinates  of  the  } joint  B in  the  two  planes,  positive  values 
having  been  assigned  to  the  radicals  between  D and  E. 

Became  of  the  singularities  at  $ — ±i  k and  at  i = ± 1,  the 
function  *1  given  by  Eq.  >; 4 1 must  be  evaluated  with  reference  to 
a path  >:  as  shown  in  Fig.  2d":<  which  has  semi-circular  Indentations 
at  each  of  these  points,  real  values  being  assigned  tc  O akag  BE. 
As  can  be  readily  verified.  Eq.  (4  - transforms  the  upper  half  cf 
the  /-plane  into  the  region  in  the  li-plane  shown  in  Fig.  2b.  In  par- 
ticular, it  will  lie  noticed  that 
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Its  the  limiting  ease.  A*  = 0.  and  the  last  formula  gives  / = see.fi. 
The  complex  potential  w can  he  expressed  in  terms  of  ( by  the 
method  of  sources  and  sinks: 

?r  = — 1 — - ^bv!  ?— / ' — In  > £ — 5 >i_| 


From  this  relation-dup  and  from  she  definition  of  L 
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The  integrals  lt  and  are  evaluated  in  the  Appendix.  If  various 
values  are  assigned  for  / and  the  corresponding  values  of  3.  «. 
and  b can  be  calculated  from  Kip.  fa),  fG;:.  and  (7.*.  The  results 
are  presented  graphically  in  Foe.  3. 

A solution  for  the  limiting  case  of  an  infinitely  deep  bucket 
can  be  obtained  without  recourse  to  complicated  analysis.  In  this 
instance.  the  velocity  k zero  at  the  bottom  of  the  cup  and  the 
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pressure  t>  correspondingly  pV  - 2.  The  moment  uni  equation  (with 
reference  to  Fig.  2a  t for  the  ./•-<!  i reel  ion  is  thus 

'2p\'rd  -p  2pF /-<!  cos  ft  - pV j~a 


a = '2  r/(1  -f-  cos  f3)  ($) 

which  is  a limiting  form  of  He;  Mit.  It  follows  that  the  maximum 
\aluc  of  u d is  4 for  this  special  case,  as  was  also  noted  in  the 
preceding  paper  in  which  this  limit  was  considered  ns  the  180° 
V-shaped  deflector.  Another  limiting  case  is  that  of  the  simple  flat 
plate  which  is  approached  as  h tends  toward  zero. 

In  experiments  conducted  hy  Mr.  Dclfaven.  observations  were 
made  of  the  angle  through  which  jets  were  deflected  by  variously 
proportioned  cups  of  cylindrical  shape.  These  arc  compared  with 
the  values  calculated  for  two-dimensional  flows  in  Fig.  8.  In  spite 
of  the  fact  that  the  experiments  were  performed  for  flow  with 
axial  symmetry  whereas  the  theoretical  calculations  are  based  on 
two-dimensional  flow,  the  agreement  between  the  results  obtained 
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from  experiments  and  those  calculated  from  theory,  equal  area 
ivitk*v  le-jug  .a-')] mod.  excellent. 

Professor  If.  T.  'Mtihly  -of  the  mathematics  depart) men t suggested 
methods  for  the  evaluat  ion  of  the  integrals  t ..countered.  Mr.  T.  T. 
Sii»  familiarized  the  author  with  the  5 * X *0 1 7*  * » a Did  5 foe  Editors 
assisted  extensive]*  in  the  preparalioai  of  the  manuseript. 


Appexdix 


The  integral  /,  can  he  devonriMw**!  into  five  eomporeiit  inte- 
grals : 
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In  which  .1 . /i.  and  C are  elementary  integrals,  and  can  be  taken 
from  Sables,  and 


» = 2(1+/)  f - --  2(  i +/  ‘ Ktkj 

\ ■ I — »-  5 1;.  1 — /:*«  * 1 


A’,  being  the  complete  elliptic  integral  of  the  first  kind.  Is  givers  in 
tables  of  elliptic  functions  for  various  values  of  k Thins  only  the 
integral  K remains  to  be  evaluated.  If  the  following  substitutions 
are  made. 
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is  ail  elliptic  integral  of  the  third  kind,  rnfortunatcly.  ihis  internal 
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»ia>  noi  uccn  taouiaieu.  so  inai  one  is  oongeu  to  express  it  m 
♦crnis  of  oi her  functions. 

If  one  writes  t = sn  (a,  A)  and  x = sn  (y,  A*),  in  which  sn  (a,  k) 
is  an  elliptic  fund  ion  of  Jacobi  defined  implicitly  by 
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in  which  sn  a is  written  ‘or  sn  (a,  A).  If  a new  quantity  c is  de- 
fined 1)Y 


so  that 
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Thef a-iK'net km  of  Jacobi.  Since  the  function  O { >j.  c)  is  an  even 
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rend  the  evaluat Losi  ©t  It  is  completed ; 
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:«  wlsick  the  phis  or  minus  sign  is  used  depending  on  whether  / is 
greater  than  i or  less  than  — 1.  Similarly,  the  evaluation  of  the 
integral  /.„  for  both  negative  and  positive  values  of  the  parameter 
/.  ran  be  amiKiplidied.  The  results  are: 
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with  tin1  same  convention  for  signs. 

For  numerical  calculations,  use  has  been  made  of  the  tables  of 
Jacobian  elliptic  functions  by  Milm-Thomson  [23 1.  An  account  of 
the  theory  of  Jacobian  elliptic  functions  can  be  found  in  [24]. 
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III.  rXSYMMETRH'ALi.V  PLACED  SEMI-INFINITE  PLATE 

by 

SlMON  lX(’K*  and  Cl.AKK  Df.1I.WKN* 

Deflection  of  a jot  by  a plate  is  a problem  in  hydraulic  engineer- 
ing for  which  a theoretical  solution  can  bo  obtained  only  in  very 
particular  instances.  The  method  of  calculation  of  the  force  exerted 
by  a jet  impinging  perpendicularly  on  a plate  of  infinite  dimensions 
is  shown  in  most  textbooks  of  fluid  mechanics.  Mathematical  treat- 
ment of  similar  problems  for  somewhat  more  complex  circumstances, 
however,  is  not  possible,  and  acceptable  approximations  must  be 


Fig.  1.  Custom  ary  Concoction  of  Jkt  Division  hy  Pki.ton-Whfki,  Bremers. 

made  unless  experimental  data  are  available.  Thus,  in  the  con- 
ventional design  of  a Pel  ton  wheel  it  is  expedient,  for  want  of  more 
accurate  information,  to  assume  that  the  portion  of  the  jet  deflected 
by  the  bucket  is  in  simple  proportion  to  the  penetration  of  the 
bucket,  and  that  the  remainder  of  the  jet  continues  in  a straight 
line,  as  shown  in  Fig.  1. 

A theoretical  analysis  of  the  two-dimensional  counterpart,  of 
this  problem  should  shed  some  light  on  the  general  problem  and 
possibly  be  of  some  value  to  turbine  designers.  Although  there  is 
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I niversity  of  Iowa.  Iowa  City. 

n 1 


. i fi.m* n*>\  of  J i : i>  — III 


oo 


no  mathematical  iust ificat ion  for  the  application  of  these  theoretica1 
results  to  actual  three-dimensional  Hows,  this  a])proaeh  seemed 
promising  in  view  of  the  fortuitous  correspondence  between  results 
of  theon  for  two-dimensional  Hows  and  those  of  experiments  for 
three-dimensional  flows  which  has  been  noted  previously  in  the 
cases  of  contraction  coefficients  for  cone  valves,  orifices,  and  mani- 
folds. 

The  conditions  of  flow  in  a Pelton  wheel  indicated  in.  Fig.  1 
can  he  approximated  if  the  flow  is  considered  to  he  two-dimen- 
sional and  the  buckets  to  be  simple  Hat  plates.  Considered  in  this 
analysis  is  the  general  case  of  a two-dimensional  jet  impinging  on 
a semi-infinite  flat  plate,  as  indicated  .schematically  in  Fig.  2a  in 
which 

a is  the  thickness  of  the  initial  jet. 

(I  is  the  perpendicular  penetration  of  the  plate  into  the  jet, 
y is  the  angle  which  the  plate  makes  with  the  horizontal, 
b is  the  thickness  of  the  guided  portion  of  the  deflected  jet, 
ft  is  the  ultimate  angle,  measured  counter-clockwise,  which 
the  free  portion  of  the  jet  makes  with  the  horizonal,  and 

V„  is  the  velocity  of  the  jet  at  infinity. 

The  ratio  of  the  division  of  flow  h /a  and  the  ultimate  angle  of 
the  jet  ft  are  to  he  defined  in  terms  of  the  penetration  d at  various 
inclinations  of  the  plate. 

From  the  energy  principle  it  is  seen  that  the  velocities  at  A. 
B,  and  E (each  of  which  is  at  an  infinite  distance  from  I))  must  all 
be  equal  for  irrotational  flow.  From  application  of  the  momentum 
equation  in  the  direction  parallel  to  the  plate,  the  relationship 

b (a  — b ) cos  (ft  -j-  y)  — a eos  y = 0 ( 1 ) 

is  obtained. 

From  Fig.  2a  a geometrical  relationship  between  b,  d,  and  ft 
can  be  found : 

d = a —J>dij  — a—  (a  — b)  cos ft  + k (2) 

in  which  the  contour  integration  is  performed  along  I)EA,  and 

A-  = fd  y + f dy 

PE  EA 

As  shown  in  the  following  calculation,  it  is  convenient  to  calculate 
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determined  from  the  momentum  relationship.  If  the  value  of  J;  for 
this  part ieular  ease  ean  he  ealeulated  in  terms  of  a.  b.  y.  and  [3. 
then  Kq.  (2)  yields  the  correspondin'!  value  of  d.  Thus  the  prob- 
lem is  reduced  to  the  determination  of  k. 

From  Fig.  2a.  dij  = sin  0 <]s.  6 being  (he  variable  angle  of  in- 
clination of  the  jet  and  s being  measured  alomr  the  contour.  Then, 
aiong  i)I\  and  EA.  Y{,d:<  = — d </>.  in  which  4>  is  the  velocity  poten- 
tial, and  one  has 


k = — 


~ sin# 
A “jr 
1)  E * 0 


d <t> 


- f 

EA 


sin# 


d <t> 


W/ 


Geometrically,  k is  the  distance  between  the  horizontal  line  through 
D and  the  horizontal  asymptote  to  the  streamline  EA  after  EA  is 
shifted  a distance  u~b  in  a direction  normal  to  the  inclined  asymp- 
totes of  I) E and  EA.  In  region  A this  shifting  results  in  the  lower 
streamline  being  moved  upwards  a vertical  distance  equal  to 
( a — b ) cos  (3.  The  integration  depends  upon  the  expression  oi  <f>  in 
terms  of  6.  and  this  is  accomplished  by  the  use  of  conformal  trans- 
formations. 

The  logarithmic  hodograph  of  the  How  pattern,  denoted  as  the 
fi-plane  ami  shown  in  Fig.  2b,  is  obtained  from  the  relationship, 

o=-ln  g-  + ie  (4) 

* 0 

in  which  q is  the  magnitude  of  the  variable  velocity  vector  and  0 
the  angle  that  it  makes  with  the  horizontal.  The  interior  of  the 
semi-infinite  strip  of  the  fi-planc  is  mapped  onto  the  upper  half 
of  the  /-plane  (Fig.  2c)  by  means  of  a Schwarz-Christoffel  trans- 
formation. the  connecting  relationship  being 

o = — x cos-1  y + i 7 (5) 

The  m-plane  is  shown  in  Fig.  2d.  The  transformation  equation 
between  the  tr-plnnc  and  the  /-plane  is  obtained  directly  from  the 
complex  potentials  for  sources  and  sinks,  point  A being  considered 
as  a source  ol'  strength  2 uYn  and  points  B and  E a*  sinks  of 
strengths  2 bY„  and  2 (<t~  b)Y,„  respectively: 

«’  = | b hi (/  — 1 ) -f-  ( n — b)  hi  (/  — e ) — a in ( / — /)  | -{-M  (0) 

in  which  M is  a complex  constant.  After  introducing  the  relation- 
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ship  dir  = d<p  + id 0 and  noting  that  along  a streamline  d 6 = 0, 
and  hence  dir  = d<p.  one  can  write  tin-  contour  integral  as 


1 , 
k -■  — j ~p—  d w 


If  dir  is  replaced  by  the  corresponding  value  from  Eq.  (6),  the 
expression  for  1;  becomes 


k = 


(7_t 

t-p 


dt 


(7) 


Also,  since  q = V,,  along  the  free  streamlines,  if  is  simply  i9  for 
all  points  on  DK  and  I) K from  Eq.  (4). 

Introduction  of  this  value  into  E(p  (fD  and  solution  for  t yield 
the  relationship, 

t = see  (y  — Q)  (8) 


which  is  valid  for  all  points  on  the  free  streamlines.  Tf  now  the 
values  of  t in  Eq.  (7)  are  expressed  in  terms  of  0 from  Eq.  (8), 
the  limits  of  integration  for  9 (y— irandO)  taken  from  the  Q-plnne, 
and  the  resulting  expression  integrated,  the  value  obtained  is 


k = -- } 6 sin  y In 


2 b 


a(l  + cos7)  — 26 


C-.  + rt  sin  7 cox  7 ’m 


2[«(l  + eos7)  — 26] 


6(1+  cos  7 ) 


— (^n  — 6^  cos  7 sin  (d+7^  ln(eos7  + cot^  sin  7^ 


(9) 


in  which 


sin  (£  + 7)  = 


\/  a ( 1 — cos  7 ) [ a ( 1 + cos  7 ) — 2 6 ] 


a 


Discussion  on  Kksui.ts 


Equation  (9)  has  been  evaluated  for  various  divisions  of  flow 
and  for  plate  angles  equal  to  20°.  4.V,  (>0°.  !)()  . 120°.  and  br>0°. 
The  corresponding  values  of  plate  penetration  and  deflection  angle 
are  presented  in  graphical  form,  with  b/a  and  /?  plotted  against 
d/a  (Figs.  3 and  4). 

It  was  noted  with  considerable  concern  at  first  that  for  very 
small  values  of  b/a,  the  penetrations  of  the  plate  d/a  were  small 
negative  values.  Since  a thorough  check  of  the  calculations  revealed 
no  errors,  further  confirmation  was  sought  by  comparison  with 


d/a 


Fi<i.  4.  Di:fi.i:ctio.\'  A.ngi.k  as  a Function-  ok  Pi.atk  Pi:ni-:tration. 
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computed.  identical  negative  penetrations  were  obtained  from  this 
indo pendent  cheek. 

The  physical  anomaly  of  a finite  division  of  flow  by  a plate 
lying  outside  t lie  path  of  the  jet  assumes  an  immediate  signiticm.ee, 
however,  if  it  is  considered  that  the  negative  of  the  complex  func- 
tion ir  also  satisfies  the  Laplacian  equations.  The  physical  inter- 
pretation of  the  negative  ir  function  corresponds  to  a reversal  of 


Fi<;.  5.  Variations  ok  b/a,  ft t and  <lm in/a  with  7 

the  direction  of  flow.  The  resulting  physical  phenomenon  in  the 
case  under  consideration  is  then  that  of  a free  jet  joining  a guided 
jet  near  the  end  of  the  guiding  plate,  the  combined  flow  continu- 
ing as  a horizontal  free  jet.  As  can  be  seen  from  the  several  graphs, 
there  is  a limit  to  the  so-called  “negative  penetration”  determined 
by  the  requirement  of  the  continuity  and  momentum  relationships. 
The  variation  of  with  7 is  shown  in  Fig.  5. 

The  behavior  uf  b/a  and  ji  in  the  region  of  negative  penetration 
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is  indicated  in  Fig.  u,  which  is  drawn  to  a greatly  enlarged  scale. 
It  can  he  clearly  seen  that  as  b/a  increases  from  0 to  0.02.  d/a 
decreases  first  to  its  minimum  value  of  —0.00243.  and  then  increases 
to  0.  From  this  point  on.  the  analysis  can  he  considered  to  repre- 
sent the  problem  posed  at  the  outset  oi'  this  investigation.  The 
physical  explanation  of  this  phenomenon  is  that  an  infinitesimal 
penetration  * causes  r.  finite  initial  angle  of  deflection  fi ; and  a 


Fig.  6.  Df.fi-kction  Anc.i.k  axi>  Division  of  I’i.ow  for  Nic.ativk 
Fknf.tr ATI () N s (y  90°). 


finite  initial  division  of  flow  b/a.  Figure  5 contains  a plot  of 
values  of  fi-,  and  b/a  as  functions  of  the  plate  angle  y. 

In  the  region  of  positive  penetration  there  is,  for  a given  y, 
an  asymntotic  limit  to  the  deflection  angle  fi  and  the  corresponding 
division  of  flow  as  the  penetration  is  increased  to  infinity  (Figs. 
3 and  4),  these  limits  having  beer*  determined  for  the  various  plate 
angles  by  means  of  the  momentum  relationship. 

An  attempt  was  made  to  verify  experimentally  the  results  ob- 
tained by  theoretical  analysis.  A jet  of  water  with  a nominal 
diameter  of  l.f>  inches  was  penetrated  by  a plate,  and  the  division 
of  flow  and  angle  of  deflection  were  measured.  Since  the  momentum 
relationship  is  the  determining  factor,  ihe  area  ratio  of  the  divided 
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jet  \v;is  substituUni  for  the*  b a ratio  of  t lit"  two-dimensional  case. 

Although  experimental  results  verity  the  general  trends  of 
the  analysis  ( Figs.  '■)  and  4 A there  are.  nevertheless,  discrepancies 
of  as  much  as  .“>()' from  the  predicted  values.  Part  of  those  can  lie 
attributed  to  t In*  crudities  of  the  experimental  arrangement,  al- 
though it  is  recognized  that  a significant  source  of  error  may  well 
lie  in  the  inherent  differences  between  two-  and  three-dimensional 
flows.  Tints  the  fortuitous  correspondence  between  two-dimensional 
theory  and  its  three-dimeus’onal  counterpart  which  has  been  en- 
countered in  other  investigations  does  not  materialize  in  this 
instance. 

However,  it  is  obvious  from  analysis  and  experiment  that  neither 
of  the  two  common  assumptions  mentioned  in  the  introduction  to 
this  paper  is  valid.  That  is,  in  considering  the  flow  into  Pelton- 
wheel  buckets,  (a)  the  division  of  flow  b/a  is  not  in  simple  pro- 
portion to  t ho  penetration  d/a,  and  (hi  the  unguided  portion  of 
the  jet  does  not  continue  in  the  same  straight  line,  but  is  deflected 
through  a considerable  angle.  Thus,  as  shown  in  Pig.  1,  the  portion 
of  the  jet  not  caught  by  bucket  „1  is  generally  assumed  to  he  caught 
by  bucket  B,  whereas,  due  to  the  deflection  of  the  jet,  it  may  actual- 
ly miss  the  wheel  entirely. 

In  view  of  the  uncertainties  of  the  applicability  of  the  simpli- 
fied theory  to  the  complex  shapes  of  Pelton-whee!  buckets,  and  in 
view  of  the  lack  of  complete  experimental  confirmation  throughout 
the  full  range  of  the  analysis,  no  positive  recommendations  can  be 
made  until  more  representative  analyses  and  complete  experiments 
are  available.  Nevertheless,  despite  the  laek  of  general  quantitative 
results  for  circular  jets  striking  buckets  on  Pelton  wheels,  the  man- 
ner in  which  a jet  is  divided  and  deflected  has  been  indicated 
qualitatively  as  a guide  for  design. 


MANIFOLD  EFFLUX 
by 

Arthur  Tocii*  and  Kobkrt  ML  Moorman** 
Introduction 

In  navigation  locks  and  turbine  scroll  eases,  water  is  often  sup- 
plied to  the  lock  chamber  or  to  the  turbine  through  a manifold 
system.  The  resultant  flows  have  been  the  subject  of  many  model 
studies  for  the  solving  of  specific  design  problems  or  the  obtaining 
of  general  information  regarding  the  flow  phenomena  involved. 
Also,  mathematical  treatment  of  manifold  flow  was  made  by 
Mc.Nown  and  Hsu  [12]  in  which  certain  phases  of  the  flow  were 
evaluated. 

A manifold  can  be  represented  in  its  simplest  form  by  an  open- 
ing in  the  wall  of  a straight  conduit.  The  lateral  can  thus  be  elimi- 
nated or  a lip  of  any  length  can  be  added.  This  representation  cor- 
responds to  some  forms  of  constructions  actually  used  [2fi]  and 
was  the  basis  for  the  analysis  by  McNown  and  Hsu.  A different 
type  of  manifold  is  also  used,  in  which  the  conduit  changes  in  cross 
section  either  in  steps  at  each  lateral  or  continuously  over  most  of 
its  length  [27].  The  analysis  in  this  paper  pertains  to  this  step- 
ped form  for  the  particular  case  of  the  infinitely  long  lateral.  Cor- 
relation of  the  results  of  the  mathematical  treatment  with  experi- 
mentally obtained  values  for  these  manifolds  is  also  discussed. 


Mathkmatical  Analysis 

A view  of  the  two-dimensional  manifold  is  shown  in  the  2-plane 
of  Fig.  1.  Flow  with  velocity  U/;  occurs  in  the  conduit  of  width  b. 
Part  of  the  flow  enters  the  infinitely  long  lateral  through  an  open- 
ing of  width  a and  attains,  after  contraction  of  the  flow  section  to 

•Roseate!)  Associate.  Iowa  Institute  et  Hydraulic  Research,  State 
University  of  Iowa.  Iowa  City. 

•Now  Associate  Professor  of  Mechanics  and  Hydraulics,  Clemson 
Agricultural  College.  Clemson.  South  Carolina. 
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Fig.  1.  Transformation  Pi.anrs. 


Cc<t  at  infinity,  a velocity  l The  remainder  of  the  flow  continues 
in  a conduit  of  width  c and  attains  an  ultimate  velocity  Tv  The 
difference  between  the  initial  and  final  widths  of  the  conduit  is 
denoted  by  d. 

Evidently,  the  evaluation  of  a coefficient  of  contraction  Cc 
and  of  a velocity  ratio  1Y  F/{  as  functions  of  the  boundary 
geometry  is  the  primary  objective  of  the  present  investigation.  The 
boundary  geometry  is  expressed  by  the  ratios  a/b  and  d/b,  so  that 
the  desired  functions  aie 
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Cc  = fx  {alb.  d!b) 

and 

Vc/Vb  = / 2 [a/b.  d/h) 

Tn  order  to  determine  these  relationships,  the  complex  function 
v must  be  defined  os  a function  of  z;  the  manifold  shown  in  the 
'-plane  is  thus  related  to  a simpler  shape,  the  infinite  rectangle  of 
the  tc-plane.  Two  intermediate  planes  are  needed  to  obtain  this 
indirect  mathematical  relationship.  The  first  of  these  intermediate 
planes  is  the  logarithmic  hodograph  or  Q-plane  of  Fig.  1.  The 
velocities  along  the  boundaries  of  the  2-plane  are  plotted  through 
use  of  the  definition  equation 

8 = - ‘n  jr.  = — In  ^ + i&  (1) 

In  this  equation  q/V j is  the  ratio  of  any  velocity  in  the  z-ph’.ne  to 
the  final  jet  velocity  in  that  plane,  and  0 is  the  angle  formed 
between  the  velocity  vector  and  the  X-axis. 

In  the  other  intermediate  plane,  the  /-plane,  the  boundary  of 
the  flow  is  prescribed  to  coincide  with  the  real  axis.  Thus  the 
Schwa r/.-Christoffel  transformation  can  be  used,  which  here  takes 
the  form 

ft  = M cosh-1  t + N 

The  semi-infinite  rectangle  of  the  ft-plane  now  becomes  the  upper 
half  of  the  f-plane.  Three  points  can  be  arbitrarily  located  in  this 
plane.  Points  A and  E are  accordingly  placed  on  the  horizontal  axis 
at  equal  distances  from  the  origin  and  are  assigned  the  values  of 
-f-1  and  — 1,  respectively,  and  D is  placed  at  infinity,  as  shown. 
The  values  of  M and  N in  the  transform  can  be  evaluated  by  substi- 
tuting simultaneous  values  of  ft  and  t for  the  points  A,  E,  arid  D , 
in  successive  steps.  Then, 

ft  = ^ eosh-1  t - i | (2) 

The  final  step  in  relating  the  manifold  in  the  2-plane  to  the 
infinite  rectangle  in  the  re-plane  is  also  accomplished  b>  means  of 
the  Schwar/.-Phristoffel  transformation, 

= *+i*  = Wirv/Hi-Vui-i)  + N' 


w 


G6 


Frkk-Stbhami  ini:  Axai.ysks 


After  integration  of  this  equation,  the  constant  .1/'  can  bo  evalu- 
uated  hv  equating  the  differences  in  the  imaginary  parts  of  the 
two  sides  of  t tie  equation  for  points  A,  B.  and  C.  The  additive 
constant  .V'  is  similarly  determined  from  the  difference  in  the  real 
parts  of  the  quantities  at  D and  E.  Introduction  of  the  values  of 
M'  and  A*'  in  t lie  integrated  transformation  vields 


bV, 


T. 

TV  \ 


c V 


C,n 


In  (/-/)  + - g)  + -c-ln(/-H)  +C-ca  (3) 


7T  Vj 


TV 


The  distances  / and  g shown  on  the  /-plane  and  occurring  in 
Eq.  (3)  can  be  expressed  as  functions  of  the  velocity  ratios.  To 
this  end,  the  expression  for  Cl  of  Eq.  (1)  is  introduced  in  Eq.  (2) 
and  the  resulting  equation  evaluated  at  point  B.  so  that 


IV  a. 
Vi*  ' 


and  at  point  so  that 


Once  the  transformations  are  completed,  the  opening  a can  be 
related  hj  the  laimal  efflux  by  integration  across  the  flow: 


a = j>  d y + Cc  a 

A E 


The  integral  can  be  evaluated  after  proper  substitutions  have  been 
made.  If  s is  defined  as  the  distance  along  the  free  streamline 
measured  from  A in  the  2-plane,  the  relationship, 


(hj  — sin  6 ds 

is  obtained. 

Because  ^ = 0 along  the  streamline  AE,  as  indicated  in  the 
tf-plane  (Fig.  I),  cw/cs  = ctp/cs  between  A and  E,  By  definition, 
d<f>/cs  = — V; ; and  dg  can  now  be  expressed  as 


dy  = 


sin  d , , 
d<p 

1 3 


From  Eqs.  (1)  and  (2), 


sin  e = + 1 ) 

along  AE.  If  the  last  two  expressions  are  substituted  into  the 
equation  for  a, 


I 


t 
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d-ccn  = fE  ~ y/^(f  + 1)  d 4> 

in  which  d<j>  can  he  expressed  in  terms  of  I by  Eq.  (21.  The  'result- 
ing integral  can  be  written  as 


a — Ce  a = 


6__.  Vh 

t \'2  Vj 


d t 


c_Vc  [ 

\/2  Vj  J+x 


Vt±} 

t—n 


dt 


Cc  a f yj+_l 
7 r v'2  J+l  f 4-  1 


df 


(4) 


in  whieh  the  limits  are  evident  from  examination  of  the  /-]>lano. 

The  evaluation  of  this  integral  is  comparatively  straight-for- 
ward, and  this  form  was  chosen  for  that  reason.  Other  substitu- 
tions than  the  above,  or  use  of  a hodograph  (£-plane),  result  in 
integrals  of  more  complex  form.  That  the  resultant  equations  will 
he  identical  to  the  one  obtained  from  Eq.  (4)  i>:  obvious.  The  tedi- 
ousness of  these  other  integrations  makes  mention  of  them  worth- 
while, however,  because  integration  for  another  How  model  may 
well  be  easiest  in  another  plane. 

The  three  integrals  in  Eq.  (4)  were  evaluated  separately,  and 
values  of  / and  ry  were  substituted  to  give  the  final  expression, 


a —Cco  = - [*  b (1/t,  4-  lVanh-1^  — 

7T  L M j-  / I j 


The  premise  of  continuity,  implicit  in  the  above  discussion,  can 
he  stated  as 


VBb-Vc(b-d)-VjCca  = 0 


(6) 


Although  Eqs.  (5)  ami  (6)  determine  the  general,  flow  field, 
another  eqiudion  is  neees;a;ry  to  impose  a particular  stagnation 
condition  which  is  physically  applicable.  The  reference  velocities 
and  the  'dimensions  of  the  boundary  should  satisfy  a requirement 
that  the  horizontal  and  the  vr.ili^ai  velocity  components  on  the 
dividing  streamline  will  vanish  simultaneously  at  D.  On  the  Q-plar.c 
(Fig.  2).  ii’  JirD  corresponds  to  the  dividing  streamline,  v will 
vanish  before  u if  the  stagnation  point  D is  approached.  Similarly, 
if  11GD  corresponds  to  the  dividing  streamline,  u will  vanish  be- 
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fore  v.  The  condition  tliat  u and  v should  vanish  simultaneously, 
therefore  requires  F,  or  G,  to  approach  J). 

The  line  BFI),  or  HOT),  is  a streamline  issuing  from  the  source 
H.  the  fi -plane  being  viewed  as  another  physical  plane  for  a fic- 
titious flow  with  F,  or  G.  as  the  stagnation  point.  Thus  at  F.  or 
G,  dw/dQ  = 0.  and  since 


dw 

dt 


d w dti 
dQ  dt 


d-  w 


dt  2 \^tJ—  1 


in  which  t is  different  from  + 1, 

dw 

dt 


= 0 


The  last  equation  characterizes  the  point  F,  or  G.  If  t approaches 
zeio,  after  substitution  of  w from  Eq.  (3)  in  the  last  equation,  a 
condition  is  obtained  which  is  equivalent  to  the  requirement  that 
F,  or  G,  should  coincide  with  l).  This  condition  equation  is 

VjCca-Vcbg-Vcdg-Vabf  (7) 

In  [12],  this  point  was  overlooked  ; it  was  called  to  the  attention 
of  the  anthers  of  that  paper  by  A.  Craya  of  Grenoble.  France.  As 
a consequence,  their  results  include  the  cases  for  which  F,  or  G, 
is  not  at  point  />.  A statement  of  this  fact  was  made  in  |28]  but 
no  analytical  result  was  presented. 

Equations  (5),  (G),  and  (7)  can  now  be  arranged  to  obtain 

Ktpf)]  <«> 


(9) 


Manifoi.ii  Ekki.cx 


69 


in  which 


Ce 


a / b 


tanh-1 


(10) 


The  values  of  the  velocities  E/{  and  Yc  relative  to  V;  in  these  equa- 
tions must  also  satisfy  Eqs.  (5)  to  (7).  The  two  functional  rela- 
tionships mentioned  earlier  as  the  primary  objective  of  the  analysis 
are  now  replaced  by  three  parametric  equations. 

Numerical  calculation  of  Eqs.  (8)  to  (10)  yields  the  graphs 
shown  in  Figs.  0 and  4.  or  comparable  ones  if  the  coordinate  axes 
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arc  chosen  different ly.  The  graphs  indicate  more  clearly,  of  course, 
than  do  the  equations,  the  How  conditions  possible  under  the  as- 
sumptions made.  In  Fig.  '•)  the  ratio  d/b  is  shown  to  be  negative 
as  well  as  positive.  The  negative  sign  denotes  mereh  that  the 
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downst roam  conduit  width  e is  greater  than  b.  As  n/b  approaches 
zero,  a limiting  value  of  for  CV  is  obtained,  as  shown 

on  Fig.  4.  The  ratio  of  the  conduit  velocities  is  the  other  parameter 
determining  larger  values  indicating  less  efflux  through  the 
lateral  and,  therefore,  greater  contract  ion. 

Experimental  Verification 

As  mentioned  in  t ho  introduction,  an  analysis  1 i 2 1 was  prev- 
iously made  of  the  manifold  with  an  opening  in  the  side  of  the 
conduit.  Comparison,  of  the  results  of  that  investigation  with  those 
of  experiments  with  an  orifice  in  one  side  of  a lock  conduit  [26] 
is  made  in  Fig.  f>.  The  curve,  in  each  of  the  two  graphs,  is  the 
mathematical  relationship  and  the  points  plotted  are  the  experi- 
mental values.  The  agreement  between  the  trends  of  experiment 
and  hydrodynamic  theory  is  remarkable,  because  the  calculations 
are  based  on  two-dimensional  flow,  whereas  the  actual  flow  was 
definitely  tbree-dimeusionai. 

Although  this  agreement  is  heartening,  it  must  be  viewed  as 
fo  .nitous.  The  attempt  at  verifying  experimentally  the  model  set 
forth  in  this  paper  was  completely  unsuccessful.  In  the  laboratory 
model  used  to  verify  the  theory,  the  coefficient  of  contraction 
could  not  he  measured,  because  a reverse  eddy  formed  at  the 
jet  surface.  The  attempt  was  made  therefore,  to  correlate  the 
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velocity  ratios,  with  the  goometi y of  the  dividing  flow.  This  proved 
to  lie  impossible,  not  much  because  the  dividing  stream  surface 
was  hard  io  determine  by  means  of  the  dye  injection  used,  but 
because  the  flow  pattern  showed  considerable  divergence  from  two- 
dimensional  conditions. 

Agreement  between  calculation  and  experiment  depends  on  the 
physical  applicability  of  the  assumptions  made  in  the  mathematical 
derivation.  Values  shown  in  Fig.  f>  were  obtained  with  free  efflux 
from  a sharp-edged  port.  The  water  jet  was  thus  surrounded  by 
air.  a condition  which  comes  close  to  satisfying  the  assumption  of 
a velocity  discontinuity  involved  in  the  free-streamline  theory.  The 
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infinite  lateral  was.  of  necessity,  a branch  conduit,  and  a reverse 
eddy  occurred  in  the  lateral  besidr  the  jot.  Although  some  experi- 
liicntal  r vidciico  ( 12]  exists  which  indicates  that  tor  high-velocity 
submerged  jets  the  eddy  does  not  invalidate  the  assumption,  for 
the  low-velocity  jet  in  this  experiment  the  assumption  was  hardly 
approximated.  The  tendency  of  the  lateral  to  draw  a dispropor- 
tionate share  of  low-velocity  fluid  from  ihe  conduit,  as  observed 
by  Barton  and  Escobar  |2S],  may  also  he  of  importance  in  this 
case. 


(Joxcuusion 

Of  die  several  possible  types  of  manifold  efflux,  two  were  con- 
sidered in  this  paper.  For  these,  the  use  of  the  free-streamliue 
theory  of  hydrodynamics  provided  a good  approximation  to  efflux 
from  a port  in  the  eond  .it  hut  failed  to  describe  the  flow  into  a 
lateral. 

The  disagreement  between  the  limited  experimental  results 
and  the  relationships  obtained  from  the  mathematical  analysis 
presented  in  this  paper  results  from  the  lack  of  applicability  to 
this  real  flow  of  the  assumptions  implicit  in  ihe  theory.  That  such 
lack  need  not  nullify  the  use  of  the  theory,  however,  is  indicated 
by  the  results  obtained  for  port  efflux.  Although  other  conditions 
of  real  flow  may  also  he  compatible  with  the  free-streamline  theory, 
the  range  of  applicability  cannot  be  determined  from  the  limited 
evidence  at  hand. 


HEAD  LOSSES  IN  MITEL  BENDS 


by 

Hanky  if.  Ambrose* 


The  sudden  deflection  of  flow  produced  by  a miter  bend  in  a 
pipe  line  causes  separation  of  the  fluid  from  the  solid  boundary 
ar  the  inner  .juncture  of  the  bend.  The  resulting  contraction  and 
subsequent  expansion  of  the  flow  downstream  from  the  bend  pro- 
duce a dissipation  of  energy  which  is  primarily  attributable  to 
eddy  diffusion.  The  head  loss  can  thus  he  approximated  by  the  con- 
ventional < ({nation  for  an  abrupt  expansion.  In  dimensionless  form 
this  is 


_h_L_ 


1--3N'2 

■Cc  V 


(1) 


in  which  /t,-.  is  the  loss  in  head,  V02/2g  is  toe  velocity  head  in  the 
pipe  with  fully  expanded  flow,  and  Cc  is  the  ratio  of  the  contracted 


area  to  the  area  of  the  pipe — a type  of  coefficient  of  contraction. 
A schematic  representation  of  the  flow  is  shown  in  Fig.  1. 

Determination  of  the  contraction  coefficient,  therefore,  suffices 
for  an  estimate  of  the  head  loss.  It  is  {lie  purpose  of  this  study  to 
determine  to  what  degree  the  contraction  coefficient  of  the  miter 
bend  in  a pipe  can  he  approximated  by  an  analysis  of  its  two- 
dimensional  counterpart 

The  two-dimensional  miter  bend  (Fig.  2a)  consists  of  three 

•Now  Assistant  Professor,  Department  of  Civil  Engineering,  L'ni- 
versitj  of  Tennessee.  Knoxville. 
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semi-infinite  planes  which  serve  as  solid  boundaries  to  the  flow, 
the  remainder  of  the  boundary  being  a curved  surface  of  separa- 
tion upon  which  a constant  pressure  is  presumed  to  act.  The  planes 
BA  and  BC  form  the  outer  boundary,  intersecting  with  the  deflec- 
tion angle  a at  the  point.  B (chosen  as  the  origin  for  z = x-\-iy). 
The  inner  solid  boundary  DA  is  parallel  to  BA  and  separated  from 
it  by  a distance  a.  The  free  surface  DC  is  a curve  and  becomes 


Lossks  in  Mii'  r Hi  nds 


parallel  to  lit'  asymptotically.  I lie  ultimate  thickness  of  the  jet 
being  b.  The  point  of  separation  1)  is  defined  by  the  abscissa  c and 
the  ordinate  a.  The  dashed  line  extending  downstream  from  D 
represents  the  position  of  the  downstream  portion  of  the  inner  wall. 

Since  the  pressure  is  assumed  to  he  constant  along  the  free 
surface,  it  follows  that  the  velocity  along  the  surface  varies  in 
direction  only.  Inasmuch  as  either  the  magnitude  or  the  direction 
of  the  velocities  along  all  boundaries  is  constant,  the  necessary 
requirements  for  application  of  the  free-streamline  theory  are  ful- 
filled. 


The  problem  can  be  considered  as  the  determination  of  the 
coefficient  of  contraction  t\-  = b/a  in  terms  of  the  angle  « and 
the  coordinates  of  the  point  of  separation  1) — that  is.  the  definition 
of  the  function  f in  the  equation 


b 

a 


(2) 


Mathematically,  however,  it  is  more  convenient  to  solve  for  c for 
gi\  en  values  of  the  deflection  angle  and  contraction  coefficient. 

The  solution  parallels  that  of  the  paper  by  Messrs.  Siao  and 
Hubbard  since  the  present  problem  is,  in  effect,  an  extension  of 
the  problem  treated  therein.  The  flow  pattern  is  first  transformed 
to  the  ordinary  hodograph  by  the  basic  relationship 


1 dw  _ _ u . v 
f "FJ'dF  “ ~V~j  lV) 


€ 


-ie 


Since  the  velocity  along  the  free  streamline  has  the  constant  magni- 
tude Vi,  DC  appears  as  a circular  are  in  the  £-plane.  The  solid 
boundaries  become  radial  lines  with  their  terminal  points  located 
as  shown  in  Fig.  2b. 

If  the  angle  a is  set  equal  to  v/n,  n being  an  integer  greater 
than  unity,  the  flow  ’pattern  can  be  simplified  further  by  the  auxi- 
liary transformation 

r=(-On  = e-”£"  (3) 

In  the  £'-p!ane.  the  entire  solid  boundary  falls  on  the  rial  axis 
(Fig.  2e ) . 

By  means  of  the  method  of  images,  the  given  flow  pattern  in 
the  t'-plane  is  obtained  directly  fr  an  the  locating  of  point  sources 
of  strength  2itV„  at  .1  and  .1'  (c'.p  - jo  7>|'M  and  a sink  of  double 
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strength  at  B.  Although  the  complete  pattern  is  symmetrica!  about 
t he  real  axis,  only  the  lower  half  of  the  £'-plane  is  relevant.  The 
eompiex  potential  is  then  written  immediately: 

?r  = 4>  + i ■ - In  ( r + 1 ) — — - In  ( — f b a ) n | 

7 r 7 r 

— — In  | r-  ( a/b )"  ] = --  - j 2 In  (c1™  f"  + 1 ) 

— In  ( e*>n  ?'-(b  n ) " 1 - In  | e',r"  fB—  ( a fb  ) ’*  ] [ (4) 


The  complex  coordinate  of  point  D in  the  2-plane  is 


cJria  = J‘rlz  — J> 

CD  CD 


1 (he  _ /»  h r 2_  _ 

V;  r ~K  ]■  7r L e iwn in + 1 


i 

ei,ni '"-(k/a)n 


F— (a/b) 


i«_ 


neiwn  (n~'  (I 


• tt  Jc  L r 

A * ~|  f 

fn  - c~ '*n(b/a)n  tn—e-i™(a/b)HJ 


h’  + c-i* 


(5) 


Equation  (5)  consists  of  the  sum  of  three  integrals  of  the  type 


['  zn~2(1z 

J o r~Kn 

in  which  K is  a complex  constant.  Since  n is  a positive  integer  the 
integration  can  he  completed  by  substituting  for  each  of  these 
terms  a polynomial  expansion  of  ?i  terms, 


r-8rff 
r‘-K  ‘ 


1 ^ p-Hr*/n 

nK  2 "(T TKT-c7*771"  d^‘K^ 


whence 


r czsji  = j_  v'1.-,.,.;.  ln  (.izi'ciL) 

J0  i"  ~K"  nK  £<*  V A'  / 


If  the  above  relationship  is  substituted  into  Eq.  (5),  simplified, 
and  separated  into  real  and  imaginary  terms,  the  value  of  c/b  is 
obtained : 
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For  « ^ 7r/2  1 ho  supplementary  angle  a'  is  defined  by 


, * * 
a — 7T  — - = — 

n m 


from  which 


m 


If  the  indicated  substitution  is  made  in  Eq.  (6),  the  value  of  c/b 
corresponding  to  an  angle  greater  than  tt/2  is  obtained  through 
aieps  differing  only  in  detail  from  the  integration  procedure  out- 
lined above: 


/2 r + 1 — m \.  ro  0 /2r-M  \“] 

= V — — V ln  L2  -2cos  r>  J 

r — 0 
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in  which 
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can  be  obtained  as  functions  of  b a.  The  results  are  plotted  in 
Fig.  3,  which  is  the  graphical  representation  of  Eq.  (2).  In  par- 
ticular, the  functional  relationship  between  c/a  and  b/a  for  a ~ 7r/2 
can  be  shown  to  be  identical  with  that  obtained  by  MeNown  and 
Hsu  f 12 1 for  the  special  case  of  total  lateral  efflux  from  a two- 
dimensional  conduit  with  a semi-infinite  barrier. 

In  Fig.  3,  for  eaeli  value  of  b/a  must  equal  zero  for 
c/a  = cot  a.  The  limit  of  unity  is  approached  by  b/a  for  all  angles 
between  0 and  tt  as  c/a  becomes  very  large  negatively.  With  a = 0, 
b/a  is  equal  to  unity  for  all  real  values  of  c/a,  whereas  for  a 
approaching  tt.  b/a  approaches  zero  for  all  real  values  of  c/a. 

Superposed  on  Fig.  3 is  the  locus  of  b/a  for  the  especially  in- 
teresting ease  in  which  point  D lies  on  the  bisector  of  the  supple- 
mentary angle  t*' ; that  is,  for  c/a  = — cota'/2.  This  case  defines 


Pig.  4.  Hkai>  Loss  in  Mm:i;  Bonus  or  V.xr.iocs  Anglos. 


so 


F rk  I.  - Sr  sr.  a M i . i n x.  Anai.  v m : s 


t h**  equivalent  contraction  coefficient  for  the  usual  design  of  miter 
bends. 

Finally.  Fig.  4 shows  a plot  of  head  loss  in  miter  bends  for 
various  angles  based  on  the  substitution  of  h/n  from  Fig.  3 for  C, 
in  Eq.  (1).  Included  in  the  figure  are  typical  experimental  results 
[29.  30,  31  j.  Values  for  the  Wcishach  curve  [32]  were  calculated 
from  tl.o  empirically  determined  equation 

yTj2g  0.9457  sin2  " + 2.047  sin4  ~ (8) 

The  agreement  of  experimental  data  with  the  curve  calculated 
by  Eq.  (1)  and  the  five-stream  line  theory  is  good  for  small  to 
moderate  deflection  angles  but  shows  increasingly  poorer  corre- 
lation for  larger  '.ingles  This,  however,  is  iu  be  expected  in  view 
of  the  increasing  importance  of  the*  secondary  flow  pattern  and 
the  resulting  secondary  losses  with  increasing  bend  angle. 
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